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The following changes to the subject report should be noted: 

Page 25, last line in third^^full paragraph: 

Change "l2 percent" to read 23 percent. 

Page 28, line 3: Change "lOOO" to read "1190." 

Page 31^ equation (l9) should read: 

w = 0 .225 JO- U„itC^m 
100 

Page 32, the first equation for M should read: 


M = 0 .225 -1- U_m 
100 ° 




Page 32, equation (2l) should read: 


M = 0.225 


200 (1 - cos 0]y[^ 


Page 32, equation (22) should read: 


w = 0 .225 


HUpmA 

200(1 - cos 01^) 


Page 38, center of page, the computation of M should be 
written: 


M = = 0.225 

A 


HU m . , 

WA ' Ah — 6.13 lb /hr, sq ft 

200(1 - cos 0 b) 

JUL 8 1948 I 


Page 39 , last line, the computation of q should read: 


q = 14.8 X 1.633 X 35 - 6.13 X 35 - 0.832 X 14.8 X 0.8(2.05)^ 
= 1018 Btu/hr, sq ft 


Page 47^ table 

VI, the last 

. three 

columns of the table should 

read . 

as follows : 




0 

0 

0 


<48 

<1.78 

<2.06 


144 

19.2 

22.3 


384 

16.9 

19.6 


720 

12.7 

14.7 

Page 51^ table 

X, column 8 

(flush 

windshield) should read: 



0 




460 




830 



1190 
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A METHOD FOR CALCULATING THE HEAT REQUIRED 

FOR WiroSHTEID THERI-IAL ICE PREVENTION 

BASED ON EXTENSIVE FLIGHT TESTS IN 

NATURAL ICING CONDITIONS 

By Alun R. Jones, George H. Holdaway, 
and Charles P. Steininetz 


SUMMARY 


An equation is presented for calculating the heat flow required 
from the surface of an inteiT>ally heated windshield in oi-der to 
prevent the formation of ice accretions during flight in specified 
icing conditions. To ascertain the validity of the equation, 
comparison is made "between calculated valu.es of the heat required 
and measui’ed values obtained for test windshields in actual flights 
in icing conditions. » 

■The test windshields were internally heated and provided data 
e.ppli cable to two common types of windshield configurations j namely 
the V-type and the type installed flush with the fuselage contours. 
These windshields were installed on a twin-engine cargo airplane and 
the icing flights were conducted over a 3.arge area of the 
United States during the winters of 1945-46 and 1946-47. In 
addition to the internally heated windshield investigation, some 
test data were obtained for a windshield ice -prevention system in 
which heated air was discharged into the ■windshield boundary layer. 

The general conclusions resulting from, this investigation are 
as follows: 

1. The amo’ont of heat required for the prevention of ice 
accretions on both flush- and V-type windshields during flight in 
specified icing conditions can be ca,lcula,ted with a degree of 
accuracy suitable for design purposes. 

2 . A heat flow of 2000 to 25 OO Btu per hour per square foot 
is required for complete and continuous protection of a V-type 
windshield in flight at speeds ^P to 3 OO miles per hour in a 
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moderate cumnlous icing condition T For the. some degree of protection 
and the same speed range, a.yalue of 1000 Btu per hour per sq^uare 
foot suffices in a moderate stratus icing condition. 

3. A heat supply of 1000 Btu per hovir per square foot is 
adequate for a fliish windshield located well aft of the fuselage 
stagnation region, at speeds un to 300 miles per hour, for flight 
in both stratus and moderate cumulous icing conditions. 

4. Ihe external air discharge system of windshield thermal 
ice prevention is thermally inefficient and requires a heat supply 
approximately 20 times that required for an internal system having 
the same performance. 


INTRODUCTION 

For several years the NACA has engaged in a broad research 
program on the problem of the prevention of ice formations on 
airplanes. Particular attention has been given to the utilization 
(in various thermal’ ice— prevention systems) of the available waste 
heat in the engine exhaust gases. 

One part of this icing research program has been concerned with 
the investigation of thermal means of windshield ice— prevention. 

The first satisfactory solution developed was the double— panel— type 
system described in reference 1 and tested on the Lockheed 12— A, 
Consolidated B-24, and C\irtiss-Nright C-46 airplanes (references 1, • 
2, and 3) • ' , ■ 

The tests of reference! resulted in the tentative specifica- 
tion 6f a heat— flow requirement, of 1000 Bttx per square foot per hour 
through the windshield outer svirf ace. .. This value- was based on data 
obtained for a V-type windshield at flight speeds up ,to I50 miles 
per hour. The flight investigation of reference 2 was also conducted 
with a windshield of this configuration, and the specification 
appeared adequate. The heated— air flush windshield installation in 
the C— 46 airplane (reference 3) provided satisfactory protection, 
but did not serve as a check on the validity of the 1000 Btu per 
hour per square foot specification for all types ,ol windshields^ 
because the windshield configuration had appreciably different icing 
characteristics than the V— types- previously tested. 

Although these initial researches provided some information on 
windshield heat requirements, the results were empirical in nature 
and could not serve as a fundamental basis for the prediction of the 


MCA TI'" No. 1434 


3 


heating requirements for windshield configurations and flight 
conditions different than those investigated. Accordingly, a 
fundamental windshield icing research was undertaken. This research 
consisted of an analytical study of the heating requirements for 
internally heated windshields, and subsequent flight tests to 
determine the applicability of the design equations resulting from 
the analysis. In the flight tests the heat flow from -the surfaces 
of various special test windshields of different conf iguratlons 
could be closely controlled and measured. Measurements were also 
made of the windshield surface temperatures, boundary-layer 
profiles, flight conditions, and meteorological factors during 
icing in order to provide data on all the basic factors which could 
be expected to affect the heating requirement. 

In addition to the research directed toward the establishment 
of the heat-flow requirement in the case where the heating was 
supplied inteimally, a secondary investigation conducted at ohe 
same time was concerned with the practicability of ice prevention 
by the means of discharging a Jet of heated air into the windshield 
boundary layer. This device was initially installed in the 
C-U6 airplane (reference 3) as a means of augmenting the internal, 
double— panel windshield system with available primary air from the 
exhaust-gas heat exchangers. The test resultvS with this initial 
installation were suff icieiitly promising to warrant further 
investigation . 

An analytical approach to the external— discharge type of 
windshield thermal ice— prevention system was attempted, but the 
unknowns involved, such as the mixing of the heated Jet with the 
boundary— layer air, precluded a reasonable prediction of the ^action 
of this system. Actual tests of external-discharge systems in 
icing conditions were considered necessary. It was hoped that 
sufficient data concerning the influence* of the various pertinent 
factors, such as Jet flow rate, size, and temperature in 
icing conditions, could be obtained to provide an indication of 
the practicability of the system and possibly form the basis for 
empirical design equations. Controlled heated-air external-discharge 
systems, therefore, were investigated for both a V— type windshield, 
and a windshield which was flush with the fuselage contours. 

The flight tests were conducted in clear air and in natural 
icing conditions with a C— 46 airplane. In addition to the windshield 
research the flight tests included the determination of wing and 
propellei' heating requirements and a study of the meteorological 
factors conducive to icing. To obtain test data in natural icing 
conditions, for all parts of the research program the C-ho airplane 
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was operated Tdj the Ames Aeronautical Laboratory, Moffett field, 
Calif., during the wirters 1945—46 and 1946—47. • ' 

The appreciation of the MCA is extended to United Air Lines, 
Inc., the United States Weather Bureau, and to the Air Material 
Command of the Army Air Forces for aid and cooperation in the 
research. In particular, the services of or James Murray, 

Avmj Air Forces, and 'Captain Carl M. Christenson and First Officer 
Lyle ¥. Reynolds, United Air Lines, who served as pilots of the 
research airplane, wore a vali.iable aid to the conduct of the 
investigation. 


SYMBOLS 

The following nomenclature is used throughout this report: 
a radius of water droplet, centimeters 


A 

Ap 

b 


c 


P 



surface area, square feet 
projected area, square feet 

pi'ojected height of windshield or half of a width of 
ribbon, centimeters 

specific heat of air, Btu per pound, degree Fahrenheit 
specific heat of water, Btu per pound, degree Fahrenheit 
radius of a sphere, feet 


B 


e 

8 

h 


J 

k 


a significant dimension of windshield used in determining 
- Reynolds number, feet 

water vapor press\ire, inches of mercury 

acceleration due to gravity, feet per second, second 

convective heat— transfer coefficient through the wind- 
shield boundary layer, Btu per hour, square foot, degree 
Fahrenheit 

mechanical equivalent of heat, foot-pounds per Btu 

thermal conductivity, Btu per hour, square foot, 
degree Feihrerheit per foot 
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c; 




m 

M 

P 

Pr 


a 

Q. 


r 


R 



""■'is, 


dimensionless quantity obtained from reference 4 and 


defined in 


the sjrmhols of this report 



length of a windshield panel or a flat plate in the 
direction of local air flow, feet 

latent heat of e v-aporati on at surface temperature, Btu 
per pound of water 

liq!’!! water content of the air, grams per cubic meter 

weight rate of water impingement per unit area, pounds 
per hour, aqmre foot 

■barometric preseiure, i'-ches of mercury 

Prandtls num'ber '' ■— — ^ > ’•ondimensional 
\ k /■ 

unit rate of heat flow, 3tu per hour, acuare foot 

rate of hea,t supply, Btu per hour 

recovery factor equal to for tur’bul&nt flow 

/ U7B A 

Ee;'r:olda number ( nondimens ional 

Vug y 

distance from the region of air stagnation, feet 
temperature, degrees Fahrenheit 

kinetic temperature rise of air, degrees Fahrenheit 

klne^tic temperature rise of the water droplets, degrees 
Fahrenheit 


u 


U 


w 


temperatiiro, degx-ees Fahrenheit absolute 

local velocity insid.e boundary layer, feet per second 

velocity, feet pei- second 

weight rate of water irapingeraent, pounds per hour 
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W 

X 

X 

y 

a 


r 

e 

4 
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weight of heated air cuppliod, pounds per hour 

distance along windshield in the direction of air flow, 
inches 


evaporation factor, defined as 1 + 

distance nom’al to the windshield surface, inches 

angle between the plane of a windshield and a plane 
normal to the streamlines around the forebody of the 
windshield, degrees 

' specific weight of air, pounds per cubic foot or grams 
pel- cubic centimerter 


0»622 Lg f ®3~®o'\ 
cpPo 


specific weight of the water droplets, grams per cubic 
centimenter 


onissivlty of windshield panel 

efficiency of water impingement, percent 

half the central angle of the total area of impingement 
on a spherical surface, degrees 


viscosity of air, gram seconds per square centimeter 
dimensionless quantity obtained from reference 4 and 


defined in. the symbols of this report asi 


^ iSy^Upb ^ 


Subscripts 

c calculated 

E experimental 

1 local conditions Just outside the boundary layer 

o reference to ambient or free— stream air conditions 

s reference to windshield external surface conditions 


av 


average conditions 
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ATTALYSIS 

Ice formations on aircraft occur duing flight through formations 
of sLipercooled droplets of vater. Upon contact with an xmheated 
surface, the water freezes and ice accumulates at a rate deperiding 
upon the efficiency with which the surface intercepts water droplets 
from the free air stream. With a heated surface the majority of 
the water evaporated is evaporated at the surface, although a few 
small drops may he evaporated in the boundary layer. The maximum 
rate of evaporation occurs when the heated surface is fuJLly wetted, 
and it is for this condition that the follcwing analysis applies. 

The unit heat flow from the outer surface of a windshield 
during flight in icing conditions can he considered as the sum of 
four individual heat losses, or 

q. = 1 2 “Is 


■where 

heat loss due to forced convection 
q2 heat loss due to evaporation of the iEipinging water 
heat loss due to warming of the impinging water 
heat loss due to radiation to the surrounding atmosphei-e 
Each of these individual heat flows will he analyzed. 


Heat Loss Due to Convection 

The equation for the heat loss due to convection, including 
the effect of kinetic heating, is written 


q = h(tg— tQ-Atjj.' ) 

CL 


( 2 ) 


where 


Atk, = 


2gJc^ 


rUi^ 

SgJCj, 


( 3 ) 
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Equation (3) is obtained from references 5 or 6 which also state 
that the recovery factor r equals Pr^/^ for turbulent flow. 


The first term of equation (3) represents the adiabatic 
heating or cooling of the air Just outside the boundary layer 
caused by the change in velocity from IT^ to Ui. The second 

term of the equation represents the heating of the air which occurs 
as the velocity is reduced from Ui to zero. Equation (3) can be 
approximated by the simpler form 


At, = 0.832 ( k ) 

100 y 


Several equations for calculating the coefficient of convective 
heat transfer h from a flat plate have been developed, and a 
comparison of these equations has been made in reference ?• For a 
turbulent boundary layer, the general form of the coefficient was 
presented as 


/ A /c-N 

h = 0.0296 X 3600 ^ O) 

j^a/=j,o.2 

Following the presentation of equation (5) iu reference various 
specific forms of the equation are derived hy expressing the 
properties of air as a function of the average temperature of the 
air in the "boundary layer. The coefficient at any distance a from 
the leading edge of a flat plate, for the region of turbulent flow, 
is given as 


where 


T 

av 

Ui 


h = 0.5X 



( 6 ) 


average temperatiure in the boundary layer 

velocity outside the boundary layer at distance s from 
the leading edge 


For a turbulent boundis^y layer extending from s-o to 
s=Z on a flat plate, the 'average coefficient is presented as 


h = 0.64 T I 


0.8 


(7) 
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Heat Loss Due to Evaporation 

As mentioned previously, the maximiijn rate of evaporation 
occurs when tho surface is fully wetted. Since this condition 
requires the maximum heat supply to maintain a given surface 
temperature., it is chosen for design. In order to establish the 
minimum heat requirement to prevent freezing of the water on the 
windshield, it is assumed that sufficient heat will he supplied to 
maintain tho temperature of the windshield surface and the water 
thereon at 32° F, Since water runhack from the heated— windshield 
area is of little consequence, no effort is made to evaporate all 
the water intercepted by the windshield; such a requirement would 
impose an exorbitant heating load in some cases. 


Using the equations of reference 5 a,nd the symbols of this 
report the rate of heat loss due to evaporation is 


qa 


0.622 h Lg 


3 — 


- I 




( 8 ) 


Heat Loss Due to Warming the Impinging Water 


Reference 5 gives the basic equation for the dissipation of 
heat to the yater that is intercepted by the windshield, but a 
kinetic heating term is added in this report since the loss of 
energy by the water droplets assvimes an appreciable value at 
airspeeds of 300 miles per hour or greater. 



where 



0.198 

2gJcp^^ \ looy 


(9) 


and M is the average weight rate of water impingement per. -un.it 
area of windshield and is given by 


M = 


0.225 


100 


U^m 


A 


( 10 ) 


Since the streojnlines and drop trajectories for various types 
of vindshields are not now available to calculate the efficiency of 
vater impingement data for spheres or flat ribbons can bo 
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utilized to give approximate impingemert values. Curves presenting 
efficiencies of watei' impingement for cylinders, spheres, and 
ribbons are given in reference 4 together with a discussion of 
their use. 


Heat Loss Due To Radiation 

The heat loss due ,to radiation from the windshield surface to 
the surroundin.g atmosphere is usually of small magnitude. It can 
be calcu].ated from the Stefan— Boltzman equation: 



VlOO/ j 


( 11 ) 


Summation of Heat Looses 

Combi’^ing equations (l) to ( 11 ) , the complete equation for the 
dissipation of heat from a windshield surface in conditions of 
icing may be written 


1 = 


h Tt.-t ^-0.832 



0.622 hL^( 

^p ^o 




r / 


T. 


.Y 


100 y \iooy 




(12) 

For flight speeds under approximately 200 miles per horn-,, 
equation (12) may be simplified by neglecting the heat loss due to 
radiation and the kinetic heating of the impinging water. Eegroupin.g 
the terms to segregate the primary heat losses, which are due to 
convective heat transfer and evaporation, from the secondary effects 
of heating the impinging water and kinetic heating of the air, equa- 
tion (12) becomes 

/u 

• q = hX(t3-to) + M(tg-to) - 0.832 


( 13 ) 
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vhere 


0 , 622 I >3 f ®s“^o \ 


(14) 


A similar grov>.ping was used ir, reference,g 5 8. 


DESCEIPTIOW OF APPAEATUS 

The researoh reported herein was corducted in flight using a 
C— 46 cargo airplane (reference 3) modified to incorporate an NACA 
thermal ice— prevention syst,en which permitted continvious operation 
in natural icing conditions. The airplane is shown in figure 1 as 
equipped for the winter of 1945—46, and in figure 2 as operated in 
the winter of 1946—47. 

Meteorological equipment was installed to measure the free— 
stream air temperature and the water drop size, droxv-size 
distrihution, and liquid water content in icing clouds . A shielded 
thermocouple connected to a millivoltmeter was used to measure the 
air temperature. Eotating cylinders of different diameters, 
mounted coaxially on a single shaft, were used to deteraine the 
liquid water content, drop size, and drop— size distrihutlons. The 
maximum drop size was determined from the area of wa.ter impingement 
on a single, nonrotating cylinder. A. detailed description of the 
meteorological instruments, their use, and tjnpical meteorological 
data obtained are presented in reference 9. 


Internally Heated Windshields 

Three general types of aircraft windshields were tested to 
provide fundamonta,! heat— transfer datai: flat— plate, flush, and 
V— type windshields. Bie first t\70 types were tested during the 
191+5—46 operations and the third during the 1946—47 opeivitions. 

Each windshield— panel test section wa,s heated by electrical power, 
and the installations arc shown in figures 3, 4, and 5« 

Flat— pl ate windshield ,- This installation was intended to supply 
heat— transfer data for tv^o-dimensional flow over a flat plate Inclined 
at differe^nt angles to the air stream. It was considered that the 
data would be applicable to the design of flat— pane windshields such 
as those common in fighter and transport aircraft. The test panel 
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VC..-J located on tho right aide of the fuselage of the airplane as shewn 
in figure 3. The panel was mounted in such a manner that tho angle 
hetween ite surface and the relative wind direction could be varied. 
Provision was made to fix tho panel in one of three positions: 30 °, 

45°, and 60 °, measured from the tangent to the fuselage at the hinge 
line. 

Hie over-all surface dimensions of the flat plate were I3 by 
18 inches, and of this area I3 by I6.5 inches, or 1.49 sauare feet, 
were uniforEfLy hec,ted. Hie heating elements were thin metal ribbons, 
and verb capable of dissipating 2500 Btu per hovr per square foot of 
panel surface when connected in a 24— volt circuit. For the base of 
the panel, a 1/2— inch sheet of fabric— laminated plastic was used, 

•The heating ribbojis wore cemented to this base and then a l/64— inch— 
thich shoot of laminated pla.stic was cemented on top to present a 
smooth, electrically insulated nuri’ace ■ to the air stream. Centrally 
located in the base material was a 4.5— i~'chr-square heat meter 
(calibrated theriiiopile , reference 10), 1/64— inch— thick, to measure 
heat flow away from the test SEirface. Electrical power was supplied 
from an alternator through a 120/24— volt transformer, and the 
voltage impressed or. tho teest section was controlled by means of a 
carbon— pile rheostat . 

Temperatures on the surface of the flat plate were measured at 
10 points by iron— con stantan thormocouples which were 0.002 inch 
thick. Tlie thermocouples were recossod in grooves 0. 003 inch deep 
in the l/64_inch— plastic outer sheet, and their junctions located as 
shown- by figure 6. The outer surface was thinly painted to provide 
a smooth surface and to protect the thermocouples. Additional 
themiocouples wore 3.ocated on the back surface of the plate and o.t 
the heat meter. 

All temperatures were recorded by a self— balancing automm.tic— 
recording potentiometer with an estimated over— all accuracy of ±3° F. 
This degree of accuracy was made possible by the provision of several 
calibration thermocouples, with temperatures Icnotm to ±1*^ F, as a 
check on the recording potentiometer. The potential of the heat 
meter was determi'^.ed manually with an additional potentiometer. 3y 
measuring the temnerattire at the heat meter, an accuracy of 
measua’enent of heat lo-sses to +10 percent was possible. 

For the determination of the boundary— layer— velocity profile 
in clear air, conventional pressure rakes were utilized and the 
pressures were recorded by photographing a laanoEieter board. 


MCA ITT No, 1434 


13 


Iqe formations on the panel were photographed vith a 
16— millimeter cinematic camera which, was remotely controlled. The 
camera was located within the fuselage about 3 feet for^^ard of the 
flat— plate panel, and was fitted with a 90^ angle lens which 
protruded into the blister showr' in figure 3* 

Flush windshield ,— ‘Ihis windshield was installed flush in the 
copilot’s panel which,. at that time (1945—46), \ras flush with the 
fuselage contours. (See figs. 3 4.) 'Ihe test surface was 

almost sq^uare, though slightly curved, and c.overed Sin area of 
1.16 square feet. Constiniction of the panel consisted of three • 
sheets of glass bonded together with sheets of plastic. Heating; ' 
elements of Tory fine resistance wires were imbedded between the . 
outer glass layer and the plastic inner layer. These., wires were 
0.0012 inch in diameter ai>d spaced 0.6 inch apart. The panel was a 
commercial product and was designed for a maximum heat dissipation 
of 1000 Btu per hour per square foot. Pow'ex* was supplied by direct— 
current generators and the impressed voltage was varied with a 
rheostat. THieimiocouples 0.002— inch— thick were cemented to the inter 
and outer surfaces of the test area and covered with clear spar • 
varnish. The locations of the thermocouple junctions- a2re shown- in' 
figure 6. 

Ice formations were photographed with a box camera having flash 
attachments and located in the pilot’s compartment. Temperatures 
and boundary— layer— rake press\u?es were recorded by the same instru- 
ments used for the flat— plate windshield. 

Y-,tyT)e- windshield .- For the 1Q46-47 operations the flush wind- 
shields were replaced by the Y— type wii'^dshield shown, in figure 5* 

The windshield center post was s.loped at an angle of 57^ (measured 
from the horizontal during cruising flight) and the included angle 
of the V was 86^ measured in the horizontal plane. The fabrication 
of the pilot’s and copilot’s panels was identical and consisted of 
three sheets of glass bonded together with interposed sheets of 
plastic (fig. 7). Heating of the panels was attained tlirough a. 
transparent, electrically conductive surface coated on the outer 
sheet of glass on the .surface adjacent to the plastic Inner lay-er. 
The panels were a commercial product and were capable o?? dissipating 
2000 Btu per hour per square foot. Althoiigh the panels were 
trapezoidal in shape, the heated area was rectangular with the dus 
bars applied the full length of the rectangle along the top and 
bottom edges to insure uniformity of heat distribution (fig. 7)* 
heated area of each panel was l4.5 ty 25 inches or 2.52 feet. 

Power was supplied by an alternator and the voltage was controlled 
with a variable trans.formor . Thermocouples 0.002— inch— thick were . 
cemented to the inner and outer surfaces of both windshields and 
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covered with clear spar vanish. Locations of the thermocouple 
..junctions are shovn. in figure 6. 


Photographs of ice accretions on the windshields were taken 
with a camera having a synchronized f3.ash attachment. Temperatures 
were recorded in the same manner as described for the flat— plate 
windshield. 


External Air— Heated Windshields 

Flush, windshield .-^ TJie original windshield heating system 
installed in the C— 46 airplane Included both an internal double- 
panel heating system and ar external— discharge heated— air heating 
system as discussed in references 3 original system 

was modified as shoT-m in figure 8 for the 194>-46 operations. The 
inner heating system was removed, the heated-air discharge slot for 
the outer system was decreased to three— sixteenths inch, and the 
extent of the discharge slot was reduced to a width of 1 foot at 
the bottom of the pilot's windshield, starting at the center post. 

The heated air was supplied by exhaust— gas heat exchangers, and the 
air flow rate was measured with a ^i?ich venturi meter. 'A butterfly 
valve was installed to control the quantity of heated air flowing 
from the discharge slot. 

Windshield surface temperatures were measured with thin iron— 
constantan thermocouples cemented to the glass and covered with 
several coats of spar varnish. Wine thermocouples were located on 
the outside surface and three were located -on the inside surface, 
as shown in figure 6. Thermocouples were located at the ventux'i 
meter and just inside the discharge slot to measure the heated— air 
temperatures. Tv/o temperature rakes were installed on the windshield 
oiiter stirface for the determination of the temperature profile in 
the heated— air jet in a few clear-air flights, Haese rakes consisted 
of nine thermocouples at different heights above the windshield 
surface up to 1—1/2 inches. All thermocouples were connected to a 
recording potentiomveter, and the estimated over-all accuracy was 
about ±3^ F, 

A pressure rake was installed for the determination of the 
heated-air jet— velocity profiles in clear-air flights. The rake 
consisted of l8 pressure probes mounted in the form of a triangular 
prism with its base on the windshield surface. Three of the 
, pressure probes were static tubes and the remainder total— pressure 
tubes. These pressures, along with the venturi meter pressures, 
were- recorded by photographing ah alcohol manometer bear’d. 


NACA TN No. 1^3^ 


15 


Two methods "wero employed to ohtain an indication of the 
quantity of heat flow fi'om the external Jot through the windshield 
into the airplane cahi?i. Ihe first method consisted of the calcula- 
tion of the heat flow based on the measured, tc.mperature difference 
between the windshield surfaces and the known thermal conductivity 
of the windshield materials. In addition to this calculation, a 
1)-. 5-inch-square heat meter was located on the inner surface as shown 
fcy'figure 8. The potential across the heat meter was read by an 
indicating potentiometer. Ihe estimated accuracy of the heat flow 
as determined with the heat-meter installation is approximately 
±15 percent. 


V-tvne windshield.- For the 19hG~hl operations with the^V-type 
windshr^ldT^ne 3/l^irich external di.scharge slot was located along 
the bottom edge of the pilot's windshield and another along the 
conter post between the windshields with the heated air 
over the pilot's windshield. Each slot was 1 foot in length. The 
1|— inch venturi meter was replaced by two 3— inch ventxiri meters 
equipped with thermocouples to measure the flow rate independently 
to the two separate discharge slots over the pilot's windshield. 

The air-flow rate to each slot could be. adjusted independently 
although a certain amount of interdependence existed, since both 
slots received air from the same source of supply. A oos ei ower 
was installed in the main supply duct to increase the oOtal quantity 
of heated air available to the slots. 

The fabrication and location of thermocouples installed on the 
pilot's panel were described in the preceding description Oj. o 
Y-type windshields. In addition to the thermocouples on 
shield surface, four thermocouples were installed just ^ 

center-post discharge slot, and six thermocouples 
inside the slot at the bottom of the pilot s windshield. ^ Thes 
thermocouples were located so as to measure the heated^ir 
temperatures just before the air was discharged 

self-halancing potentiometer was used to record the temperc.turc . 

The pressure differentials at the yentui’i meters were indicated 
by standard differential pressvu>e gages and recorded by a standard _ 
NACA 60-cell recording manometer. 


TEST PROCEDURE 


The clear-air flights 
Laboratory, Moffett Field, 
the Laboratory served as a 


were conducted at the Ames Aeronautical 
Calif, in the case of the icing flights, 
main base' and operations in search of 
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icing conditions were conducted over a co’isiderable area of the 
United States. The 1945—46 operations, consisting of over 40 
flights, were conducted mainly in the Pacific- IJortlwest and Ilorth 
Central States. The 1946—47 operations Incivided this same area 
and, in addition, some icing conditions were encountered in the 
Central Southern states. 

The flights in clear air were made to check or modify the 
theoretical equations for convective heat— transfer coefficients 
and to measure the degree of kinetic heating experienced hy the 
windshields. 


Icing Flights ■ , 

The tost procedures employed were substantially the samo for 
each test windshield. Therefore, the following discussion applies 
to the operations of both winters. When jcing conditions were 
anticipated, the windshield power was txirned on or, in the case 
of the heated-air jet, the valves were opened. Before entering 
the icing condition, the heat to the windshields was adjusted to a 
setting greater than the expected requirements. After a stabiliza- 
tion period of several minutes the recording instruments were started 
and the cloud was entered. After the completion of one nun, if no 
ice formed, the heat supply wouD.d be decreased and the equiipment 
allowed to stabilize for 10 minutes before the next set of data was 
taken. The heat input to each windshield was reduced in this manner 
to approximately the minimijm amount required to keep the heated 
surface free of ice accretions. For the pilot's and copilot's wind- 
shields, the value of minimuim heat required was established by 
visual observation of the start of ice accretions on the windshields. 
Since the surface of the flat— plate panel could not be observed by 
the windshield engineer, the outer— surface temperatures were 
reguilated by changing the heat input so that the lowest temperature 
would be Just above freezing. 

The pilots always made an attempt to hold the altitude and 
airspeed constant duning a specific test run, but severe turbulence 
sometimes made this impossible. The duration of a test run varied 
from a few m.inutes in the top of a cumulus cloud to severc-l hours 
in an extensive stratus layer. For a few runs in the smaller cloud 
formations, the airplane was flown in, a circle in order to remain in 
the clotid long enough to obtain a complete set of meteorological and 
heat— transfer data. Uone of this latter type of data has been 
pi-eaented unless the conditions were - fairly stable for at least 
5 minutes after an equipment— heating stabilization period of 10 
minutes. 
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In addition to the flight and meteorological conditions^ a 
continuous record *was kept of the following windshield data: heat 

supplied^ inn.er— and outer— surface temperatures, photographs of ice 
formations, and the heat losses inward as indicated hy heat meters. 


Clear— Air Tests 

Clear— air flights were conducted in level flight at various 
altitudes and airspeeds to cover the conditions encountered during 
the icing operations. The test flights were made in a direction 
away from the sun to minimize the effects of , solar radiation. 

First, flight data were taken without any power supplied to 
the windshields. By measuring the surface-temperature rise above 
free— stream temperature a basis was esteblished for comparison 
between the indicated recovery factor of equation (4) , and the 
theoretical value 'Pr^/^ for ’turbulent flow. 

A check for possible significant edge losses was made by varying 
the heat supplied to the windshield panels while flying at a constant 
airspeed and altitude. Any variance in the calculated external heat- 
transfer coefficient after correction for inward heat losses \?as an 
indication of the edge losses, since the effect of changes of 
temperature in the boundary layer on the heat— transfer coefficient 
can be considered- negligible for the temperature ranges concerned. 

For several airspeeds and altitudes, heat was supplied to the 
windshield panels while the temperatures, heat losses, and velocities 
were recorded. These data were necessary for the experimental 
determination of the .convective heat— transfer coefficient h, and 
tho establishment of effective values of’ s or I (eqiuitions (6) 
and (7)). 

During the clear-air flights after the 1945-46. operations, a 
survey of the boundary layer at the center of each test windshield 
was made with a pressure rake. The data were used to deteimine the 
boundary— layer thickness and the variance of the local velocity 
outside the boundary layer with changes in airplane velocity. Both 
pressure and temperature boundary— layer surveys were made in the 
heated— air Jet over the flush windshield installation. 
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RESULTS OF TESTS 
Internally Heated Vfindsliields 

Local air relooitles and convective heat— transf er coefficients .— 

The values of local air velocity just outside the boundary layer at 
the center of the flat— plate panel and the flush windshield are 
presented in figure 9 for various airplane speeds and altitudes. In 
figures lOj 11, and 12 the convective heat— transf er coefficients, 
based on the measured stirface temperature at the center of the flat- 
plate and the known heat flow, are presented. 'The coefficients for the 
flush panel are not presented because of difficulties experienced 
with the measurement of the surface temperature. In the case of the 
V— type windshield only a few values of the heat— transfer coefficient 
obtained in clear air are available. These data were secured 
between the icing encounters and, since the data were i-ecorded at 
different flight altitudes, they do not form a satisfactory basis 
for curve plotting in the same manner as figures 10, 11, and 12. 

This information is presented in the f ollovring table : 


1 ' 

{ 

•Jtest 

no. 

Fressure, 

altitude, 

(ft) 

True 

airspeed 

(mph) 

1 

Free— air 
temperature 

(°F) 

1 

Windshield 

surface 

temnerature 

•(Of) 

h 

Btu/hr, 
sq. ft, 

(Op)' 

1 

12,000 

178 

20 

73 

23.4 

2 

10,800 

200 

18 

65 

28.5 

3 

9,8oo 

174 . 

39 

73 

23.0 

4 

7,8oo 

175 

30 

66 

26.5 

. 5 

6,000 

170 

35 

70 

25.2 ! 

i 


Meteorolopiical and heat— transfer data during icing conditions .— 

Thbles I and II present the meteorological data for the two winter's 
operations which have been selected for discussion in this report. 
These data represent only a small portion of the meteorological 
information recorded, but have been selected as the tests which 
supply the most satisfactory combination of simultaneous meteoro- 
logical and windshield heat— transfer data. A complete tabulation 
of all the meteorological data obtained in the 1945—46 operations 
is presented in reference 9, which also includes a detailed 
discussion of the procediires employed to obtain the test results. 
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The rotating— cylinder method utilized to measure the liquid we-ter 
content and meanr-effectiTe drop— size values, presented in tables I 
and II, provided a 1— minute average of the icing conditions, since 
this was the time of exposiure of the cylinders. ^ The windshield ■ 
heatr- transfer data required about ' 2 minutes to record and the 
rotating cylinders were exposed at some time during that interval. 
The meteorological data of icing conditions 4, 5j and 6 of table I 
were obtained prior to the installatioii of the rotating cylinders, 
hence the maximum drop— size range encountered during the heat- 
transfer test interval is presented. 

The heat— transfer data for all three of the internally heated 
windshields in conditions of icing are presented in table III, All 
the heat flow values presented in the table do not represent the 
minimum requirement for the particular icing condition, since cases 
for which satisfactory meteorological and heat— transfer da.ta were 
obtained together did not always occur at the minimum value of heat 
flow. Ihe surface temperatures presented in table III are the 
arithmetical average of the recorded values for the several thermo- 
couples on the surfaces. Photographs and sketches showing the ice 
accretions on the heated windshields and sxtrrounding surfaces for 
several of the conditions encountered are presented in figure I3 

to 23. 


External Discharge Windshields 

Temperature and velocity gtradients .— The tempeiutures of the 
windshield surface at various distances from the heated-air discharge 
slot for the flush windshield installation are shown in figure 24. 

The temperature and velocity gradients in the jet for the same 
windshield installation are presented in figures 25, and 26. 

Icinn; tests .— Table IV presents the windshield (flush) surface 
temperatures moasaired during operation of the external discharge 
jet in icing conditions. An indication of the ice— prevention action 
of* the external discliarge Jet for the flush windshield installation 
can be seen in figi.ire I5. .-lie effect of the V— type windshield 
external jet (pilot's side of windshield) is shown in figure 20. The 
incomplete removal of ice accretions shown in this figure was typical 


^Mean— effective drop size for a cloud sample is defined as the 
ddometer of a water drop for which the amoun.t of water existing in 
water dropvS larger than that drop equals the amount o.f water in drops 
smaller than the drop. 
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external discharge system on the V-type 
w.cndshi8ld. The temperature limitations of the windshield panel 
precluded the use of additional heated air in an attemnt to remove 
the residual accretions. 


discussion 

Predaction of Heating Roquiroment for 
Internally Heated Windshields 

Since the main objective of this report is to establish 
fundamental design equations for the calculation of windshield heat 
requirements in specified icing conditions, equation (13) will be 
discussed in the light of the test data to establish its usefulness 
for tnis pumose. The applicability of this equation for the 
determination of heating requirements can best be evaluated by a 
com^rison of calculated and measiu-ed values of heat flow from the 
windshield surface, for the same surface temperature, and for 
specified flight and icing conditions. Before such calculations 
can be made, however, means must bo established for the determina- 
tion of two components of equation (I 3 ); ncaaely, the convective 
heat-transfer coefficient h and the weight-rate of water impingement 

following discussion considers first the 
utilization of the flight-test data to establish means for the 
evaluation of h and M with a degree of accui'acy suitable for 
esign p'urpoaea. Although the methods discussed are directly 
applicable to the test windshields of this report, an effort was 
made .o generalize the technique in order to provide a design 
proceaure applicable to other configurations and icing flight 
conditions, 

th e . gQry^^Ctlve hea+^transfer co^fCi _ pj, ^he 
deteraination of the value of the convective heat-transfer coeffi- 
cient h, the pertinent question is the suitability of equation (6) 

-or the calculation or the coefficient. A comparison of measured 
and computed values of the convective heat-transfer coefficient h 
or the outer surface of the three internally heated test wind- 
shields therefore will be made, 

^ Considering first the evaluation of h for flat-plate panels, 
ligures 2f, ^.8, and 29 present a comparison of the measured values 
of h at each panel angle and at various altitudes and airspeeds 
wi .h values calculated by equation (6) . The curves representing the 
measured values of h have been taken from figiures 10, 11, and 12. 
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Eqi:>ation (6) is used for com-parison rather than equation (7), sinco 
it has heen derived for the purpose of calculating the coefficient 
at a given point on the surface and the test data apply to the center 
of each panel. The values of Ui used in equation (6) for the 
calculation of the coefficients presented in figi;res 27 ^ 28 ^ and 29 
were the local velocities at the center of the panel as given in 
figure 9* Two values of s wore used in the calciilations si 
representing the distance from the panel hinge to the ce-nter of 
the panel, and s^ the distance from the center of the panel to 
the stagnatio;! point at the nose of the airplane fuselage. 

In the case of the V—t;^’pe wl*'^.d3hield, the calculated values 
of h from equation (6) are compared with the test values (see 
table presented in Results) in the following table: 


— 

Item 

no. 

measured ' 

1 

1 ' ' 

!h calculated 

(equation (6)) 

3l = 0.67 ft ' 

' 32 = 5.0 ft 

1 

*1 

21.3 

1 ■ ■ ' 

29.1 

19.5 

2 

2p.9 

33.2 

22.1 


20.8 

29.6 

19.8 

4 

24.0 

32.4 

21.6 

5 

22.7 

33.1 

22.1 


In the computation of h by equation (6) the value of the local 
velocity over the windshield was not known and hence the airplane 
velocity was used. 

An examination of figures 27, 28 , and 29 table just 

presented indicates that the Uwse of equation (6) provides reasona-bly 
accurate values of the convective heat— transfer coefficient at the 
center of the test windshields. In terming the accuracy of equa- 
tion (6) as ''reasonable,” consideration has been given to the fact 
that no data other than the length and location of the pa?-\el and the 
local velocity were employed and the agreement is considered reason- 
able for such an approximation. A more exact determination of the 
value of h would require the application of somewhat lengthy computa- 
tion methods to measured values of the boundary— layer profile, or 
the installation of a heating plate in the windshield surface and the 
procurement of actual test data. Tlie disadvantages of these moans 
must be weighed against the inaccviracy of the general application of 
equation (6) in determining the method that will be omplo^/'ed to 
evaluate h in a^-^y fvittire desigp. coimputations. The agreement 
shown in figures 27, 28 , and 29 between the measured values of h 
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and the values computed equation (6) is about the same^ 
irrespective of the value of s employed. In most instances • 
the two computed cum^es bracket the experimental curve ^ witli the 
lower value of s providing the larger values of h. For design 
purposes the use of the lower value of s should provide heat— • 
transfer coefficients which would probably be soi!iewhat larger 
than the true case^ but would be conservative. 

The spacing of the heating wires in the flush windshield 
(0.6 in/ between wires) resulted in considerable nonuniforraity of 
the surface temperature as can be observed in t}ie ice-removal photo- 
gi'aphs (fig. 15)- As a result the atterapts to calculate the 
convective heat— transfer coefficient for the surface based on the 
moasured heat flow and surface— temperature rise produced question- 
able and erratic values. Consequently a comparison of measured and 
computed values of h for the flush windshield has not been presented. 

In an attempt to establish the a.ccuracy of equation (I3) for 
the prediction of the heat requirement^ it is desirable to utilize 
the most accurate available values of the various components of the 
equation. Thus although the foregoing discussion indicates that 
equation (6) provides reasonably accurate values of h (which can 
be used when more precise methods are not available or are 
considered too complicated for the degree of accuracy desired)^ a 
more accurate determina,tion of h would bo desirable for use v/hile 
checking the general equation. One possible solution would be to 
utilize the test values obtained in clear air_, but this system has 
the disadvantage of inflexibility since corresponding clear— air 
data were not available for all the icing tests. The alternative 
employed was to determine the value of s in equation (6) which 
would provide the best average agreement between calculated and 
measured values of h for the test range of airspeeds and altitudes. 

Ihe determination of an average s for the case of the flat- 
plate panel at an angle of ^5^ will be used as an example of the 
procedure followed. In figure !!_, 12 test points are plotted as 
measured values of h. Iho data cori*e spending to each test point 
(h^ Tj U, and 7) were inserted in equation (6) and the value of s 
computed for each point. The arithmetical average of the 12 values 
of s was calculated to bc:i 1.63 feet. For the determination of the 
value of h at the center of the flat— plate panel when set at 45^ 
for different velocities and altitude s^ equation (6) becomes: 


h = 0.51 Tav 




0.8 


(1.63) 


o .2 


= 0.46 Tav""-^ (Ui7 )^-® 


0.3 


( 1 ?) 
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This same system was employed for the other angles of the flat- 
plate panel, and for the V— type windshield, with the following 
resi:lts: 

Panel angle = 30°j s = 0.37 foot, and h = 0.62 Tav°'® (1(5) 

Panel angle = 60^", s = 1.93 feet, and h = 0.45 Ta/'® (1?) 

Y— type v?indshield: 

o = 1.26 feet, and h = 0.4l Tav'^*® (iB) 

Evaluation of the rate of water Imp ingement.— Assignment of a 
value to the factor M of eqmtion (I3) for specific design purposes ' 
is a problem regarding which very little information- is available. 

The most recent and extensive treatiaent of water--drop trajectories 
around several generalised objects is preaenued in reference. 4. This 
investigation was undertaken primarily to provide a method, for 
estimating the rate of water impingement on an airfoil and, there- 
fore, most of the report concerns the interception of v;ater drops 
by cylinders which are assumed to represent the forv/ard section of 
the various airfoils. Some trajectory calculations, however, are 
presented, for the cases of a ribbon (flat pla.te norma], to the 
direction of air flow) and a sphere. These calculations are used 
in this report to predict the rate of vra-ter Impingement on the flat- 
plate panel (ribbon) and the flush windshield (sphere). 

Considering first the case of the flat— plate panel, the 
assumption was made that the rate of vrater impingement on the panel 
would be equal to the rate of impingement on the pi-ojectod area of 
the panel considered as one— half of a ribbon. (See fig. 30- ) The 
efficiency of vrater impingeaient r, and the wei^it rate of water 
impinge,^:ient w for various ^op sizes have been computed for the 
three panel -angles at one flight condition and ere presented in 
table V. It is of interest to note that for drop d-iamators greater 
than 3C1 microns, as the panel angle is reduced (panel becoming more 
flush with the fuselage) the impingement efficiency is increased but 
the weight rate of ws.tor interception is decreased. Since this . 
latter quantity is the item of greater significance from a heat- 
requirement standpoint, the desirability of ].ow panel angles is 
evident, , Optimum benefit of this advantageous featur’o would be 
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e:cpected to occur with flush windshield installations. The flush, 
windshield in the C— 1-6 airplane during the 1945—46 winter provided 
confirmation in this respect, since wind.shield icing was encountered 
in only 8 out of 30 actual flights in icing clouds. Figure l6 shows 
the lack of ice accretions on the flush windshield after a flight in 
a cloud which produced the accretions sho^vn on 3/8-inch rods located 
just below the pilot’s and copilot’s windshields. 

A second item of interest which can be noted in table V is the 
fact that the maximum calculated efficiency of water impingement is 
8o percent, and this value corresponds to exceedingly large drops 
for an icing condition. This result is somewhat surprising, and 
possibly in error, since it would appear that the presence of the 
fuselage forvrard of the windshield would cause a concentration of 
drops near the surface and produce impingement efficiency values of 
large magnitude, possibly greater than 100 percent. Some attempts 
were made to obtain an indication of the rate of water impingement 
on the C— 46 windshields diuring the operations, but nothing satis- 
factory resxilted. 

To obtain an estimation of M. the rate of water impingement 
On the flush windshield for specific design conditions, the 
assumption was made that the impingement would be the same as that 
on a portion of a sphere having a diameter eaual to the fuselage 
maximum diameter. (See fig. 31.) Ihe manner in which the water- 
drop trajectory calculation method for a sphere presented in 
reference 4 was app].ied to the configuration shown in figure 31 is 
discussed in detail in Appendix A of this report. Tiie calculated 
rate of water impingement on the flush test panel, based .on the 
procedure discussed in Appendix A, is presented in. table VI for 
various drop— size diameters. The same flight conditions employed 
for the flat— plate water interception calculations (table V) were 
again assuiaed to apply. 

An. interesting conclusion res\fLting from an inspection of 
table VI is that, for the conditions presented, drop diameters of 
50 microns are required in order to achieve a value of 0]y[ of 
sufficient magnitude to reach the test panel. This conclusion is 
in agreement with the observations of water impingement on the 
C— 46 flush windshield in icing conditions. The drops were observed 
to strike the windshield only during flights for which the 
meteorological data indicated the presence of moan effective drop 
size of 30 microns in diameter, or larger. The fact that the rate 
of vrater impingement on the flush panel as presented in table VI 
decreases with increased drop size is the .result of the assumption 
(Appendix A) that the water is uniformly distributed over the area 
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of iiaplngoncntj and laeans that the area of inpingemont increases at 
a more rapid rate than the weight rate of water impingement. This 
possihlo error in the va.lue of the weight rate of water impingeiuent 
is considered to he of negligible magnitude in view of the suiall 
.amoimt of water involved for the drop— size range assumed. 

Establishment of appli c.?. bility of eq.'aation (I 3 ) f or the 
prediction of windshield he at ing requireme nts.— The applicability 
of equation (I 3 ) will now be demonstrated by a comparison of calcvuLated 
and measured values of heat flow in conditions of icing for the flat- 
plate, flush, and V— type windshields, liao assumptions and method of 
calculation followed are presented in Appendix B. To illustrate the 
procedure employed, sample calculations are presented in detail for 
the flat plate and flush windshield for icing condition 6, tables I 
and III. An illustrative calculation for the V— type windshield is 
not included, since the procedure followed was identical to that 
presented for the flat— i^late panel. 

The calculated values of heat flow are compared vrith the measured 
values for the flat— plate and V— t;)ip)e installations in tables VTI and 
Till, respectively. The fact that the calculated values of q are 
lower than the mcjasurod value scin table VII and higher in table VTII 
may be attributed, in part, to the different velocity basis used for 
the determi.nation of h in oaoh case. The calculated values in 
table VII should bo more nearly correct since the velocity employed 
was the local velocity ovei- the pa.ncT, while in the case of the V— type 
windshield (table VIII) the airpl:ino speed was used since the local 
velocity over the windshield was not knoifm.. 

The few instances in which icing of the flush test panel 
occurred px'ovided very limited data for a comparison of the calculated 
■and measxired values of the hea.t flow. The one case presented as a.n 
illustration in Appendix B represents the only reliable data, and 
the calculated vaTuo was 12 percent greater than the experiUBntal 
value . 


One possible reason for the la,.rgo disagreement, in some cases, 
between the measured and calculated values of heat flow nxy be 
traceable to the assumption, used for 'all of the calcvTations, that 
all the drops were of the same size. Tlxo calculation of the heo,t 
flow from the flat— plate panel for icing condition 12, table VII, 
will bo used to illustrate the influence of drop-size distribution on 
the heat .flov/. The calculated vaTue of 390 3tu pe.r houj?, square 
foot presented in table VTI, icing condition 12, wa.s based on*?, 
uiu.foi-'m drop size of 19 microns. The r'otating-cylinder data, however, 
indicated a drop--size distribution defined as tjpe E in reference 4. 
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Ife.'ble IX presents the restilts of calculating the veight— rate of water 
impingement on the flat— plate panel, aasitming the water in the free 
stream exists in drop sizes corresponding to distribution E. On this 
basis the original weight rate of water impingement (1.91 lb per hr) 
is increased to 8.14 pounds per hour. Although this is a sizeable 
Increase in the amount of water intercepted, only q.3 (equation (9)) 
of equation (I3) is affected by the increa-se, and the revised value 
of q becomes 1035 Btu per hour, square foot, which, in this case, 
provides better agreement with the measured value. 

The assumption of uniform drop— size distribution in the calcula.- 
tion of q may also be the cause of the large disagreement between 
the measured and calculated values of q in icing conditions 1 and 
7, table VII. Based on the assumption of no drops present larger 
than the mean effective diamoter, the weight— rate of water impin-ge— 
ment in these two cases is found to be zero. The value of q in 
table VII, therefore, has been derived for these two conditions on 
the assumption of a clear panel with no evaporation occurring on 
the surface. If drops of a d.iamcter large enough to strilce the panel 
had actually been present in sufficient qi;antity to wet the panel, 
even partially, the calculated value of q would be considerably 
increased and better agreement with the measured value would result. 
Although these discussions are by no means conclusive, they do 
indicate the imiportance of drop— size distribution and maximum drop 
size. This factor should be given consideration in future design 
calciilations. 

Based upon an examination of the data in tables VTI and VTH, 
the calculation for the flt^sh windshield heat flovr in Appendix B, 
and a consideration of the various influencing factors which cause 
disagreement betvreen the calcixLated and measured values of heat flow, 
it is concluded that eqviation (I3) will provide predicted values of 
heat flow with an average accuracy of 15 percent if the entire surface 
is assumed or known to be completely iretted and the local velocity 
over the windshield panel is employed in the calculations. The 
accuracy of the equation is greatly dependent upon the accuracy of 
the values of convective heat— transfer coefficient employed. The 
determination of experimental values of this quantity for the 
particular windshield configuration proposed would be of considerable 
aid in the accurate prediction of the required heat flow, although 
in the absence of such data equations (6) and (7) will supply a 
reasonable approximation. The effect on the accuracy of equation (13) 
produced by an error in the estimation or calculation of the weight 
rate of water impingement is not large, provided sufficient water is 
intercepted by the windshield to completely wet the surface. Rates 
of water impingement less than the value required to wet the surface 
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affect the valvis of qs (equation (&)), since evaporation over the 
entire surface is not realized. For rates of impingement above the ■ 
quantity required for completely wetting the surface, the va3.ue of 
qs remains constant and the only term of eqiiation (13) affected is 
the value of q3 (equation (9)) which is an item of secondary 
magnitude. ... 

(lenei^l conslderatlors of vindfihleld heating requ,lremen.t3 

Having shown that equation (13) is applicable for the calculation of 
windshield heating requirements, provided the various components of 
the equation can be evaluated, it is of interest to utilize the 
equation to investigate the heat requirements for various windshield 
configiu-ations in icing conditions which might be selected as design 
requirements. An indication of the meteorological conditions, 
corresponding to typical, or average severity, icing and to the most 
probable maximum icing to be expected in all weather transport 
operations was obtained from a review of reference 9. 3'or each of 
the two general cloud types (cuuulus and stratus) a liquid \rator 
content was selected which corresponded to a moderate icing condi- 
tion. In the case of the stratus cloud, the value assxmied was 
0.5 gram per' cubic meter, while the corresponding value for the 
cumiilus cloud was 1.0 gram per cubic meter. Based on the data, in 
figures 6 and 9 of reference 9j values of 15° F, free-air temperature, 
and 15 microns, drop diameter, were selected for the stratus cloud. 

The corresponding values for the cumulus cloud were 0° F and 20 
microns . 

The heating requirement for the flat— plate, flush, and Y— type 
windshields in the above icing conditions was calculated by the 
method presented in Appendix B for airplane speeds of I50 and 
300 miles per hour, and the restilts are nresented in table X. From 
this table the previous empirical heating requirement of 1000 Btu 
per hour per square foot of windshield surface (reference 1) is 
seen to provide adequate protection for Y— type windshields in 
moderate stratus for flight speeds as groat as 3OO miles per hour. 

In the case of the moderate cumulus cloud, however, this quantity 
of heat flow would not be adequate even at a speed of I50 miles an 
hour. In transport aircraft operations, the most critical need for 
windshield ice protection occurs during flight in the congested area 
surrounding the air terminals. Any icing condition encountered would 
be of the stratus type which would tend to reduce the heating 
requirement. For transport operation in mountainous areas, whore 
cumulus and stratocumtilus predominate and high cruising altitudes 
(hence low froe-air temperatures) are the rule, a heat require- 
ment of from 2000 to 2500 Btu pei* hoiir per square foot, for speeds 
up to 300 miles per hour, is indicated if complete and continuous 
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protection is dosirfed. Certainly complete visibility at al3. times 
is desirable for military aircraft. It is of interest to note that 
a heat supply of 1000 Btu per hour per square foot appears adequate 
for a flush windshield in all of the icing and flight conditions 
assumed in table X. 

In considering the individual contribution of each of the 
ve.riou3 meteorological factors which constitute icing conditions 
to the total heat requirement, the powerful effect of changes in 
the value of free-a,ir tempera. ture is noted. Ihis dominating effect 
is traceable to the fact that the surface evaporation is greatly 
Influenced by temperature changes, and this term qs of the general 
equation is one of the majoi- components. As an illustration of this 
effect, consider the increase in heating requirement (table X) for 
the fiat— plate windshield from 5^0 to 10^4-0 Btu per hour, square- 
foot as a result of (l) increasing the drop diameter by 5 microns, 

(2) doubling the liquid water content, and (3) decreasing the free— 
air temperature by 15° B. iJhen equ.ation (13) is utilized, it can 
be shown that the 5*^0 Bt-a per hour per square foot increase is 
composed of a 13— percent increase due to increased drop size and 
liquid water content, and 87 percent duie to the change in the free- 
air temperature, This fact leads to the conclusion that once the 
windshield heat requirement ha,s been established for a specified 
icing condition, and when the condition is assumed to be severe 
enough to completely wet the windshield surface, changes in drop 
size and liquid water content will not change the heati;ig require- 
ment apnreciably, A change in free-air temperature, however, will 
have a very noticeable effect on the heating requirement. It should 
be noted that this conclusion is at variance with wing heating 
requirements because, in the case of wings, runback of the impinging 
water is not desirable and, therefore, sufficient heat is supplied 
to evaporate all of the water wherever this is practicable. An 
indication of the amount of hoat that would be required if this 
design procedure were applied to windshields can be obtained with 
equation (I3) . Calculations with this equation, based on icing 
condition 6, table I, -irdicate that evaporation of all the water 
impinging on the flat— plate panel in this condition would require 
a hoat supply of 19,400 Btu per hour, square foot. 


External Discharge VTirdshields 

In the introduction of this report the statement was made that 
meastirements were taker, in the jet of the heated— air external- 
discharge system In the hope that a basis for the establishment of 
empirical design' equations- would result. Unfortinately, a review of 
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the velocity and temperature profiles measured in the Jet, of which 
figiires 2^ and 26 are typical examples, did not reveal any basis for 
a rational analj'-sis of the mixing of the Jet with the boundary— layer 
air, and the prediction of the resultant surface— temperature rise. 

A few items of interest^ however, were noted dtiring the invest!— 

, gation. One of these items was the large amount of heat supply 
required for ice protection by the extenial discharge method in 
comparison with internal heating of the surface. For example, 
figure 15 illustrates a flight in which the windshield area under 
the discharge Jet was maintained clear with a heat supply of 20,000 
Btu per hour, which is an approximate unit flow of 10 ,‘000 3tu per 
hour per square foot of cleared surface. In the same icing condi- 
tion, ice accretions were removed from the flush test panel with a 
heat flow of 5^5 Btu per hour per square foot. (See fig. 15 S'Pd 
table III.) Thus the external discharge systqm required a heat 
supply approximately 20 times that required for the internally -heated 
system, for the same degree of protection in the same icing condition. 

The thermal inefficiency of the external— discharge system is 
apparently the result of rapid mixing of the discharge Jet with the 
cold boundary— layer air, with a resultant rapid decrease in. the Jet 
and windshield surface temperature as shown in figure 25. The 
sujrface temperature data in this figure indicate a decrease from 
l6o^ F at the discharge slot to a value of 42° F at a distance of 
only 6 inches from the point of discharge. 

Satisfactory operation of the external— discharge system for 
the V— type windshield was not obtained because of failure of the Jet 
to flow across the entire s’jirface. Figure 20 illustrates a typical 
ice— removal attempt, A small area along the bottom of the pilot's 
windshield, and also at the center post, has been cleared by the Jet 
but the remainder of the panel is covered by an ice formation. 

Attempts to raise the surface temperature by increasing the flow 
rate and temperature of the heated air were limited by the temperature 
re'strictions of the vinyl plastic at the lower edge of the windshield 
(region of maximum temperature). The externo.1— di scharge system of 
windshield ice prevention appears to be a desirable installation only 
in those instances where (1) internal heating is not possible, 

(2) the discharged air will flow over the windshield without requiring 
the additional weight penalty of blowers, and (3) a large supply of 
heated air is available in the region of the windshield to bo 
protected. 
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CONCLUSIONS 

The following conclusions are based on the analytical studies ^ 
and test data of this report and should he applioahle ^ to windshield 
configurations and icing conditions siinile.r to those in.vestigaued. 

1. The coefficient of conwective heat transfer for the e^rtejmal 
siurface of a Y-tj^e or flush windshield- can he approxtoated with 
accviracy suitable for design purposes by the use of the. established 
eq^uations for tiurbulent flow on a flat plate. 

2. The heat requirement for ice prevention on a flush or 

Y— type airplane windshield during flight in specified icing condi- 
tions can be calcuJLated to an accuracy of 20 percent. 

3. The complete and continuous prevention of ice accretions on 
the surface , of a Y-type airplane windshield, for flight in modemte 
cumulous icing conditions at speeds up to 300 miles per hour, wj .11 
require a heat flow from the surface of from 2000 to 2500 Buu per 
hour per square foot of stirface. In the case of continuous flight 
for the same speed range in moderate stratus conditions, a heat 
flca^ of 1000 Btu per hour per square foot should prove adequate. 

h. The complete and continuous prevention of ice accretions^ 
on the surface of a flush-type airplane windshield located well aft 
of the fuselage stagnation regio-, for a speed range up to 300 
an hour in stratus and moderate cumulous conditions, can be obtained 
with a heat flow of 1000 Btu per hour per square foot of sirface. 

5. The tendency of ice to accrete on windshields which are 
installed flush with the fuselage contoirs is considerably less than 
that for T— tsrpe windshields. 

6. The external discharge sj>-3tem of windshield thermal ice 
prevention is thermlly inefficient and requires a heat sxipply 
approximately 20 times that required for an internal system having 
the same performance. 

7. Windshield installations which conform to the fuselage 
contours are more adaptable to the ^.^se of the external discharge 
system than V— tj'pe installations because the heated Jet Wj.11 flow 
naturally over the windshield surface. 


Ames Aeronautical Laboratory, 

National Advisoi-y Committee for Aeronautics, 
Moffett Field, Calif. 
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APPENDIX A 

Mejthod Used for the Determinatior of the Pate 
of Water Inplngsment for the Flush Test Panel 

Ihe cross sections of the C— 46 airplane fuselage consisted of 
two circular sections separated hy the floor line, the upper section 
heiiTg the greater. The entire flush windshield assembly was included 
between two stations which were located 23.5 ^.nd 50-5 inches, respec- 
tively, from the nose of the fuselage and are designated as stations 
A and B, respectively, (fig. 31). flnsli tost panel was located 

in-the copilot's windshield as shown by figures 3 2,nd 4, 

A sphere represents the best forra for which streamlines and 
water-drop trajectories are now hnovm (reference 4) and which would 
have, approximately the same I'ate of water impingements as the front 
of the fuselage. For calculations a sphere was selected with a 
radius equal to the maximum radius of the fuselage. The relative 
location of the flush test panel on. this sphere is shown, in figure 31. 

Water intercepted bv the snhere .— The projected area of the 
sphere is used with the impingement efficiency to determine the weight 
3?ate of water interception: 


w = 0.225 — - Uom2irC 
100 


(19) 


where 

C the radius of the sphere^ feet 

Weip;ht rate of water intercepted per unit of iirrplmement area ,— 
The weight rate of water intercepted hy the sphere divided hy the 
surface area over which the water droplets impinge gives the avera.ge 
rate of impingement over the area of impingement. 

The equation for the impingement area is: 

Ai = 2TtC^ (1 003 (20) 


where 

0^ half the central angle of the total area of impingement on a 
spherical surface 
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Therefore, the weight rate of vrater intercepted per unit of impinge- 
ment area may he written: 

r\ 27TC 

M = 0.225 Uq m 

100 2TrC^ (1 — cos %) 


which may he reduced to the form. 


M = 0.225 


■1 


100 0(1 - cos %{) 


( 21 ) 


We ight rate of water intercepted hy the flush test panel .— 

The assumption was made that, for water-drop sizes of sufficient 
diameter to cause the area of impingement to incltide the flush test 
panel, the water intercepted per unit area of the panel would he 
approximately equal to the average weight rate of water intercepted 
over the area of impingement. On the basis of this assumption, the 
follovang equation was written: 


w = 


0.225 


r| Up m A 

100 C(1 — cos %) 


( 22 ) 


where 

A the area of the flush test panel, l.l6 square feet 

Calculated values • of water impingement on the flush test panel 
for various water-drop sizes are presented in table VI. 
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APPETIDIX B ■■ 

CJilculations of the Heat Esq_uii’.ements on- the Flat— Plate 
arjd Flush Test Windshields for Icing Condition 6 

I., .Heat reguire-aient for the flat-ulato vindshield -.— 

A. The' follcw.i.-ng assumptions were made: 

1, The water intercepted by the windshield is' the same as that 
intercepted by a ribbon with an area equal to the 
; : projected area of the windshield. (See reference. 4 and 

' ■ figure 30.) 


2. The water that is intercepted by the windshield panel is 

heated to the surface temperature and then all, or in 
part, evaporated. ' ' ' 

3. The windshield surface is completely wetted. 

it-. There is a region of air stagnation at the loading edge 
of the panel. 

5. The air flow over the windshield is turbulent. 


B. 'The conditions for illustration are the same as icing condi- 
tion 6, table I : 

1. Pressure altitude: 12,700 ft 

2. Airplane true airspeed: lit-0 mph 

3. Freo-atream air temperature: —2° F 

i|. Liquid water content: ,.1.0 gram/cu meter 

p. Water— 'ircp diamotor: L4 microns 

■6, Drop— size diatrib'uti'oh : uniform 


7. Heated area of panel: 1.49 sq ft 

8, Convective heat— transfer coefficients based on test data 

taken at the center of the panel and the average value 
of s = 1.68 ft. 


9 . 
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9. Panel angle with the fuselage: 49 ° 

10 , Moan surface temperature of the panel: 36° P 


C. External conTective.. heat— transfer coefficient. Since the' 

].ocal velocity vas measured only at the center of the fLo-t— 
plate panel and .equation (Ip) (Discussion) gives the ■ 
experimental convective hoat— transf er coefficient at the 
center of the panels a conversion factor vas reeded to 
convert the coefficient at the center of the panel to an 
average value, for the entire heated area. By substituting 
2s for I in equation (7) and equating the two eqmtions, 
the average coefficient was found to he 1.1 times the 
coefficient at the center of the panel. Therefore-, by using 
equation (15) fo^r this case: 

h = 1.1 X 0.46 (t'lr)'^*^ 


T; 


av 


-t- 46 o = 477O F absolute 

2 


= 106 X 1,467 = 156 ft/sec 

7 ^ 0.0525 Ib/cu ft 

therefore • 

h =17.3 Btii/hr, sq ft, °F 
D. Bate of water impingement, equation (lO) : 


M = 0.225 


100 




At 


The percentage of water intercepted is obtained from curves 
presented in reference 4 after values of K and cp have been 
calculated from the equations of that reference: 


K = 


cp = 


2y,..^^Uo 

i8r^Uob 


= 1.21 


= 1.48 X 10' 
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where 

7^ specific weight of the water drops, 1 graia/cai cia 

7 specific weight of the air, 8,22 x 10 grams/cu cm 

a radius of the water drop, 0.0022 cm 

h projected height of windshield or half of rihhon width, 

1.5 X 30.48 X, cos 45° = 32.3 cm 

Uq airplane speed, 6,250 cm/sec 

g acceleration due to gravity, 98O cm/sec, sec 

p ahsolute viscosity of air, 1,69 x 10 gram scc/sg cm 

For these values of K and cp the impi?igemert efficiency is 
equal to 43 percent, and the equation for the rate of water 
impingement becomes : 

M = 0.225 X .43 X 205 X. X 1 = 14 lb /hr, sq ft 

149 

w = MA = 20.9 lb /hr 


E. Calculation of the factor X, equation (l4) ; 

T - 1 + Q ♦ 822Lg ( e^-ep) 

*^p ^o 

For the sample conditions 

Lg = 1074 Btu/lb at tg = 38° ? 

c = 0.240 Btu/lb, °F 

P 

= 18.55 in. Hg 

e., = 0.212 in. Hg 

Op = 0.040 in. Hg 
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tg = 36° F 

t = -2° ■R' 

k-o <= - 

'Iherefore 

X = 1 + 0.675 = 1.675 

F. Heat reauired at the outer surface of the flat— plate panel for 
ice preventj.0'4, equation (I3) : 

a = hX(tg-to)TM(t3-t^)-0.832 h 

= 17.3 X 1.675 X 38 + lA x'38 - 0.832 X IT . 3 X 0.89 (2.05)^ 


= 1580Btu/hr, sq ft 


G, Temperature of the panel surface required to evaporate all the 
vater intercepted: 

Eqiiating the equation of mass tranefer (equation (8) without 
the latent heat of vaporization) to the weight rate of water 
impingement^ the difference in vapor pressures, required to 
evaporate all the water intercepted is first determined: 

h ! ^62^ \eg-eo) =20.9 Ih/hr 


As a first approi^rimatior, h can he taken as the previously 
calculated value of 17.3 Btu/hr^ sq ft, 

then 

= 5.8 in. Eg 

e- = 5.8k in. Hg = lk8.3 mm Hg (saturation) 


therefore, 

tg = 139° F 
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For this value o; 


3g ft, 


tg, h is recalculated to be 17. S 3tu/hr, 


Therefore 


e = 6.01 in. Hg = 1^3 Hg 

s 


H. Heat required at the outer surface of the panel to evaporate 
all the water intercepted with no water runback or blovr— off . 

h =17.8 Btu/hr, sq ft, ^ 

M = li^.0 Ib/iir, sq ft 

tg = l4l° F 

X = 7.0 

q = 17.8 X 7.0 X lUi + Ik.O X l4l - 56 = 19,400 Btu/hr, sc. ft 

II, Heat reoulrement for a flush.--t.Tpe windshield .— The method 
of calculation and the assumptions made were similar to those used for 
the flat— plate windshield. Icing .condition 6, table I, is again 
selected for illustration. For the flush— type windshield the 
assumption was made that the water intercepted by the windshield is 
equal to the water intercepted by a portion of a sphere with a 
radius equal to the maximum radius of the upper half of, the fuselage. 
(See Appendix A, and fig. 3I.) 

A. External convective heat— transfer coefficient. Since the 

flush windshield test panel did not deviate from the fuselage 
contours, the distance from the stagnation regioji at the 
nose of the fuselage to the center of the panel was utilized 
in equation (6), to determine the convective heat— transfer 


coefficient. 



= k75° F, absolute 

Ui =190 ft /sec 

7 = 0.0525 Ib/cu ft 

s = 5.0 ft 
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Therefore, 

h = lh.8 3tu/hr, sq ft, °F 

B. Bate of vater impingement, eqmtion (22), Appendix A, Tne 
percentage of water intercepted was obtained from cm-ves 
prti seated in reference 1+ after the following values of L 
and cp had been calculated; 

K = 0.2k2 

Cp =7.5 X 10^" 

then 

T| 3 percent 
and 

a,, > 16° 

The central angle 0j,{ = l6® does not represent an impinge- 
ment area extensive enough to include all the area of the ^ 
flush windshield (fig. 31). yet test observations indicated 
that during icing condition 6 ice tended to form over the 
entire windshield surface. For this reason,. 63 = 27.. 5 

was substituted for in equation (22) . The weight r^-ite 

of water impingement was determined by the equation: 

T>U^ m , 

M = w = 0.225 2 = 2.3 Ib/hr, sq ft 

A 100 C (1 —cos 85) 

C. Calculation of the factor X, eqmtion (1^) 


X = 1 + - 


0 . 622 Ls /f’s-eo^_ 


^0 ^s“'^o 


tg = 33° F 
= -2° F 

L. = 1075 Btu/lb 


= 18,55 in. Hg 


1 + 0.633 = 1.633 
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Cp = 0.24 Btu/ro, °F 
©3 = 0.1878 in. Hg 
= 0.04 in. Hg 

D. Total heat required at the outer surface, equation (13) • 
q = 14.8 1,633 >; 35 + 2.3 x 35 - O.832 x l4.8 x 0.8(2.05)^ 

- 884 Btu/lu', sq ft 
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TABLE I.- ^^ETEOROLOGICAL AMD FLIGHT DATA, C- 1+6 AIRPLAME, WINTER 


Icing 
• condition 
nUDibei" 

Flight 

ntmiber 

Bate 

Time 

Preesiire 

altitude 

(ft) 

'Iru-e 

airspeed 

(raph) 

Liquid 

’(■rater 

content 

(gm/m°) 

Mean- 

effective 

drop 

diameter 
(microns ) 

.Ambient - 
air 

tempera— 

tux-e 

• (°F) 

Remarks 

1 

31. 

3/29/46 

1116 

8,000 

•160 

0.2 

10 

14 

ciEiiulus clouds 

2 

: 31 

3 / 29/46 

1628 

11,650 

180 

.5 

18 

17 

cxmrulonirilbus clouds 

3 

33 

3/31/46 j 

1147 _ 

8,300 

3.70 

.9 

10 

18 

cumulus clouds 

‘ 4 

23 

3 / 2/46 

1454 

10,350 

135 

.7 

^■11 to 20 

6 

cmulus clouds 

: 5 

23 

3/2/46 

1516 

11,500 

145 

.8 

^"20 to 30 

-.5 

— ■ ■ -- — ■ - 

cumulus clouds i 

; . 6 , 

23 

3 / 2/46 

1548 

12,700 

140 

1.0 

^18 to 44 

-2 

cumulus cloxids | 

V : 

39 

4 / 8/46 

1413 

9,000 

160 

.5 

11 

24 

stratus and strato-'] 
cumulus clouds j 

8 

39 

4 / 8/46 

1451 

8,700 

160 

.2 

13 

26 

stratus and strato— 
cvauxilus clouds ! 

9 

46 

4/25/46 

1.656 

10,300 

170 

.2 

48 

20 

aitostratus clouds ! 

largest drop size [ 
incurred 

10 

49 

4 / 29/46 

1121 

9 ,l 4 o 

165 

.7 ■ 

23 

9.5 

cumulonimbus clouds 

11 

h 9 

4 / 29/46 

1132 

8,930 

165 

. 2 , 

20 

11 

cumulonimbus clouds 

12 

49 

h/29fk6 

1155 

10,600 

160 

1.4 

I9E 

5 

cumulonimbus clouds 
largest liquid water 
content incurred 


Maximum drop diameter 
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TABLE II.- METEOROLOGICAL AMD FLIGHT DATA 
C -46 AIEPLAUE, WHTER 19k6-k'] ■ 


Icing 
condl— 
ti on 
nuniber 

Flight 

niffliber 

Date 

Time 

■■ ’ 

Pressure 

altitude 

(feet) 

True 

air 

veloc- 

ity 

(mph) 

Liquid 

water 

conteiit 

(gm/m^) 

143 an 
effect— 
tive 
drop ■ 
diameter 
(microns) 

Ambient 
air ' 
ternr- 
pera— 
tur© 
(Of) : 

Mazimimi 
drop size 
(microns) 

D 1 s — 
tri— 
hu- 
tion 

Remarks 

13 ’■ 

78 

2/2/47 

1057 

6., 780 

163 

0.17 

14 

11.2 ■ 

- 

D 

Cum'olus • 

14 i 

" 83 

2/5/47 

1243 

19,500 

150 

.22 . 

24 

-12 

- 

C 

Alto stratus 

.. ..^.4 

- 95 : 

3/10/47 

1630 

11,500 

167 

.64 

19.2 

23 

- 

A 

Cumulus 

16 i 

95 

3 /l 0/^7 

i 64 o 

11,000 

165 

.44 

16.4 ■ 

23 

- 

A 

Ciimulus . . . 

17 

100 

3/16/47 

1155 

10,840 

l 64 

.12 

10.4 

21.1 

13 

C" 

Stratus 

18 

105 

3/21/47 

1102 

5,080 

• 162 

.39 

13 

18.5 

l 4 

B 

Sdrat ocumul us 

19 ' 

105 

3/21/47 

1114 

5,000 

.152 

.34 

15 

19.3 

19 

B' 

Stratocumulus 

20 

109 

3/21/47 

1429 

5,150 

162 

• 57 

17 

21.0 

27 

E 

St rat oc ijinulus 

21 

105 . 

3/21/47 

1433 

4,980 

156 

.37 

19 ' 

21.2 

19 

D 

Stratocumulus 

22. 

105 

3/21/47 

1505 

5,190 

168 

.41 

20 

20.8 

18 

E 

St rat oc 

23 

105 

3/21/47 

1512 

5,150 

151 

.30 

22 

19-3 

23 

E ' 

Stratocimilus 

1 

f 

4 

105 

3/21/47 

1528 

4,500 

. 161 

.26 

12 

22.0 , 

16 

— 

E 

Strat oc umulus 
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TABLE III.- RESULTS OF FLIGHT TESTS OF IKTERHALLY HEATED WINDSHIELDS, C -46 AIEEJUIE 


I Surface | q 

Condition Banel I tempera- * (Btu/hrj 
nuiriber | an gle jture (°F)i /sq ft) 


WINTE R 


1Q1|5 - k6 


Flat-Plate Panel 


FD.ueh Tyne 


Surface 


8 


30 


50 


1520 


30 


39 


t 


12l^0 


30 


60 


1^5 


1230 


tempera— (Btu/hr 
jture C^F ) j /sq ft ) 


Remarks 


WINTER 1946 - 1)7 

,V^e 


37 


bondition 

number 


28 


274 'Heat supply adequate I3 


~ I iF.p.p.heafsupply. li 
2(^4 adequate ice previ-l 
! ously^f ormed remov— j! 
ing “^F.W. clear m 


Durface 1 q j 
W ind— tempera— j ( Btu/hr j 
shield ture j/sq ft) j ' 


Remarks 


Co- I 
•pilot's j 


45 


1740 


Pilot's j 27 


1030 


31 


37 


910 


45 


4o 


16.30 


45: i 36 


! 1720 


23 


.F.P.P.heat supply jj 

adeauate sli^t iceji 
I on F.V 7 . without j; 

j hea t supply j]_ 


15 


I Co- I 
i pilot ' s I 


46 


935 


Heat supply 
adequate 


Heat B'apply 
inadequate 


Heat supply 
I adequate 


280 I Heat supply adequate;! I6 


27 


565 


! Heat supply adequate 
; vioibio runback on 


P Co- 
I pilot' 3 


4 ? 


17 


! Pilot's! 4 l 


920 


.Heat supply 
I adequate 


43.0 


F.P.P. 


iHeat supply 
adequate 


33 


60 


59 


1220 


60 


-48 


640 


50 


780 


I Heat supply inade- 
i quate large j^form- 
a t i cns on d o tn 
vindshlelda 


?18 


Heat supply adequate 


18 


! Pilot’ si 38 


630’ 


19 


j! 

» I 


49 


j Heat supply not 
I quite adequate to 
I F.P.P. ShW. clear 


20 


Hlot's! 39 


Co — 

pilot ' 3 
~Uo^- 


39 


Heat supply 
i sligjitly 
' inadequate 


gc!o Heat supply 
^ { sligJitly 

i inadequate 


gjQ IHeat Bupply 
I adequate 


10 

6 c 1 

[ ! 

47 ! 

' r 

11-70 

4 ; 1 

1 26 

h- 

0 

: : M 

Heat supply adequate! 

1 — — - 

22 

'.Co- 
pilot' s 

39 ' 

6S0 

1 i 

iHeat supply 1 
1 adequate • 

11 

1 ! 

1 60 

i 

i 

1 1 

53 

1 

1190 

! 

27 

1 

0 

Heat supply not i 
quite adequate j 
ice ‘’being formed | 
and removed j 

23 ' 

1 Oo— 1 

j pilot’ s 1 

i i 

1 

1 

; 38 

! 680 ! 

1 

1 

i 

j Heat supply 
1 adequate 

(. 

12 

t ■ ■■■. ■■ 

1 60 
i 

40 

1 

1235 

1 

1 — 

t . 21 - " 

1 

i_t_ 

0 

Heat supply adequate' 
ice sliding off SP..P. r 

24 

i 

Co- 

• pilot' s 
1 

1 

i 490 

L- 

1 Heat supply! 
i subtly j 
1, AitaAeqwtg ! 


IF.P. F.— Plat— plate panel 


^F.W.— Flush— type -windshielda 
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T/lBLE IV*- RESULTS OF FLIGHT TESTS OF A EEATED-AIR 
EXTERN/i DISCH/iRGE SYSTEJ'5 FOR I'TIiroSHIELD ICE 
PROTECTION. C-46 MRPL/J^TE, l-'flNTER 1945-46 


^ICIRG CO^IDITION 4 

Heat supply slightly in- 
adequate • Windshield 
clear over three-quarter 
area of heated jet. 

Air-flow rate, 390 Ib/hr 
Nozzle air temperature, 147^ F 
Nozzle air velocity, 105 mph 
Heat supplied^, 13,200 Btu/hr 


X 

( ) 

"^s 

(°F) 

^s“^o 

(OF) 

0 

147 

141 

3 

75 

69 

6 

52 

46 

9 

41 

35 

12 

36 

30 

15 

32 

26 

18 

32 

26 

21 

31 

25 

24 

30 _ 

24 


ICING CONDITION 5 

Heat supply adequate. 
Windshield clear in the 
area of the heated jet. 

Air-flow rate, 460 Ib/hr 
Nozzle. air temperature, 187^ F 
Nozzle air velocity, 140 mph 
Heat supplied, 21,000 Btu/lir 


] 

X 

(in.) 

(OF) 

^s-^o 

(OF) 1 

. i 

0 

187 

187 

3 

96 

96 

6 

71 

71 

9 

54 

54 

12 

47 

. 47 

15 

40 

40 

18 

40 

40 

21 

39 

39 

24 

34 

34 


ICING COivTDITION 6 

Heat supply inadequate' 
for complete clearance 
of winds Meld. 


Air-flow rate, 470 Ib/hr 
Nozzle air temperature, 190^ F 
Nozzle air velocity, 150 irmh 
Heat supplied, 22,000 Btu/hi" . 


X 

(in.) 

(op) 

1 

ts~to ■ 
(op) 

0 

190 

192 

3 

83 

85 

6 

62 

64 

9 

46 

48 

12 

40 

42 

15 

33 

35 

18 

30 

32 

21 

25 

. 27 

24 

22 

24 
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^ Conditions presented in table I. 

^ Available heat in jet above ambient-air temperature 
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TABLE r?.- CONCLUDEI), RESULTS OF FLIGHT TESTS OF A HEATED- 
AIE EXTERNAL DISCHARGE SYSTEM FOR WINDSHIELD ICE 
PROTECTION C-k6 AIRPLANE, WINTER 1945-^6 


ICING COIiDITION 9 

Ho at supply adequate . 
Windshield clear. 

Air— flow rate, 470 Ih/hr 
Nozzle air temperature, 192^ F 
Nozzle air velocity, I 38 mph 
Heat supplied, 20,000 Btu/hr 


ICING COITOITION 11 

Heat supply adequate • 
Windshield clear. 

Air— flow rate, 33 0 Ih/hr 
Nozzle air temperature, 178 ° i‘ 
Nozzle air velocity, 86 mph 
Heat supplied, 13,000 Btu/hr 


ICING CONDITION 12 

Heat supply adequate . 
Windshield clear. 

Air-flow rate, 3 00 Ih/hr 
Nozzle air temperature, l82° F 
Nozzle air velocity, 86 mph 
Heat sxipplied, 13,000 Btu/hr 


I X 

1 ; 

! ■*^s 

! 

1 

{ 

! ’’^3-^0 
i (OF) 

- 1 

1 o^’l 

1 

! ^92 

i 172 1 

. ( 

1 , 

1 3 : 

: 62 

1 

42 j 

! 

1 

1 

i 55 

i' 

! 1 

! 3 ? ! 

1 1 

1 

1 

^3 • 

’ 23 

i ' 

37 

! 17 

1 15 

! — 

' - 

18 

3 ^ 

1 

14 i 

■ ' ! 

21 

1 

i 33 

♦ 

; 13 

j 24 

I 30 

; 10 


! X i 

(in.) 1 

1 - 1 

j 

( 05 .) 1 

L 

^s"^o 

m ! 

i i 

0 j 

CO 

• 168 1 

1 

f 

. 3 

70 ! 

i 

’ 60 j 

6 

i 

53 .i 

43 1 

1 ,. '! 

1 : 1 

' i 

32 ! 
1 

"■' '12 . ' 

I 38 

i ' 

28 i 

I ^5 

1 

26 

18 

i 

i 28 * 1 

j ! 

18 i 

21 

• ,,,31:.- 
1 ’ 

21 ; 

'24 

1 35 ‘ 

25 1 


' ^ 1 
1 (in.) i 

j 

<°f) 1 

^s-^o i 

(°F) ' 

1 0 ! 

1 

1 

182 ,l 

180 : 
♦ 

? 1 
; 3 1 

79 ' 

! 

77 ' 

1 i 

! 6 1 

' 57' I 

1 

55 

9 

■32 i 

30 i 

12 

1 

38 1 

36 i 

■ 15 

1 3*. 1 

32 : 

j 

18 

! 33 ' 

; f 

31 j. 

i 21 

1 30 

28 1 

1 24 

1 31 

29 1 
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TABIJ: V.- CALCULATEID VAI.IIES OF WATER INTEDRCEP'i’ED BY A FLAT-PIATE PAIOJL 
BASED ON DATJi CALCIHATED FOE EIBBOWS, EEFERFJICE k. 


-p- 

U\ 


Drop 

diainoter 

microns 

axio"^ 

2 8 

a xlO 

Panel angle 60 ^ 

Panel angle 45^ 

Panel angle 30^ 

K 

(percent) 

(Ih/hr) 

K 

. 4 

(ih/iir) 

K 

1 

4 

(iD./h^) 

].0 

5 

25 

0.097 

0 

0 

0.07 

0 

0 

0.1 

0 

0 

20 

10 

100 

.185 

1 

.64 

.227 

2 

1.04 

.3^ 

5 

1.84 

30 

15 

225 

.51 

15 

9.6 

.629 

20 

10,4 

.884 

• 32 

11.8 

50 

25 

629 

1.41 

49 

28.8 

1.73 

90 

26.1 

2.44 

60 

22.1 

100 

50 

2500 

5.65 

j 

! 46.6 

6.93 

75 

39.0 

9.8 

80 

29.5 


Pressure altitude, 10,000 f"'-. 
Tr'.ie airspeed, 1^0 nph 
AialienWiir temperatui*-?, 0® F 
Liquid water content, 1.0 ,'p7i/m.'''' 


Area of panel, 1.49 sq 
Half width of rihhon = 
Average value of cp = 2 


ft 

panel lerigth 
X 10^ 


CCS a 
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TABLE VI.- CALCULATED VALUES OF WATER EvIPINGEMENT 
ON A ELU 5 E WINDSHIELD PANEL. BASED ON DATA 
CALCULATED FOR SPHERES, REFERENCE 4 


Drop 

diameter 

(microns) 

i 

I n 

(percent) 

1 

1 

! 

» 

f 

i 

1 

1 

I 

1 rnn: — ; 

j Flush 

! teat 

I Sphere panel, 

! w M 

1 (Ih/hr) (Ih/h:^ 

1 i sq ft) 

1-. 

Flush 1 
test 
panel, 
w 

(lb /hr) 

20 

0 

■ ,.^c 

1 

0 io ■ 

I , 

0 

10 

<1 


<18 1 0 

0 

50, 1 

1 

16'^ 

54 \l .2 

8.4 

70 1 

! 8 i 

ip - I 

i 26° • 

l 44 i6.3p 

7.4 

100 

i 

; 15 

, - 1 

45° 

1 ; 

270 |4.76 



Pressure altitude ^ 10,000 ft 
True airspeed, I50 mph 
Amtient-air teinperature , 0 ° F 
cp = 1.1 X 105 

= half of the central angle 
of the total ai’ea of impingement 
on a spherical surface 
Liquid water content, 1 0 gm/m,^ 
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TABLE VII.- COMPARISON OF CALCULATED AND ^EASUEED 
HEAT FLOW FROM THE SURFACE OF THE FLAT-PLATE 
PALEL. C-46 AIEPLANE, VUNTEE 194-5-46 


-!=■ 

00 


Icin-g 

Condi- 

tion 

Number 

Panel 

iingle 

From 

Fuse- 

lage 

Mean 

Si.irface 

Temper- 

ature 

r 

Calcu- 

lated 

T| 

(50 

Calcu- 

lated 

V 

(lb /hr) 

Calcu- 

lated 

h 

(Btu/hr, 
sq ft/°F) 

Calcu- 

lated 

qc 

(Btu/hr, 
sq ft/oF) 

Experi- 

mental 

cLe 

(Btu/hr, 
sq ft/°F) 

q.c-iE 

X 100 
4 c 

1 

30 

50 

0 

0 

27.9 

890 

: 1520 

-630 

-71.0 

2 

30 

39 

9.5 

2.1 

30.4 

1210 

1240 

-30 

-2.5 

3 

30 

60 

C 

0 

30.2 

1180 

1230 

-50 

- 4.0 

4 

45 

37 

3-5 

1.15 

17.7 

730 

910 

-180 

- 24.6 

5 

45 

40 

20.5 

8.26 

18.3 

1454 

1630 

-176 ^ 

- 12.1 

6 

45 

36 

39 

19.0 

17.3 

1580 

1720 

-i 4 o 

-8.9 

7 

6 o 

59 

0 

0 

16.0 

500 

1220 

-720 

-144 

8 

60 

48 

0 

0 

16.0 

66 c 

640 

20 

•3.1 

9 

60 

53 

45 

6.44 

14.9 

1250 

1205 

45 

3.6 

10 

60 

47 

5-5 

2.7 

14.8 

1100 

1170 

-70 

- 6.4 

11 

60 

53 

2,5 

0.35 

14.9 

1300 

1190 

no 

-^.5 

12 

60 

40 

2.0 

1.91 

15.6 

890 

1235 

-345 

-38.8 
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TABLE Till.- COMPAEISOW OF CALCUL/i.TED AlJD MEASURED HEAT FLOW FROM THE SURFACE OF 
THE V-TTPE WINDSHIELD, C -46 MKPLAEE. WINTER 1946 - 47 . 


Icing 

condi- 

tion 

nuiuber 

Windshield 

Mean 

surface 

tempera- 

ture 

(opy 

Mean 

inoide 

surface 

teLiuera— 

Calcu- 

lated 

(per- 

cent) 

Calcu- 

lated 

( 14 / 

hr) 

Calcu- 
lated h* 
(BUi/hr, 

Calcti— 
la ted 

(Btu/hr, 
sq ft) 

Experi- 

mental 

(Btu/hr, 
sq ft) 

X 100 

Q.C 

Remarks 

13 

Copilot's 

45 

73 

0.1 

0.04 

34 . 5 ' 

1900 

1740 

8.4 

Heat supply 
adequate 

14 

Pilot ' s 

27 

24 

9.0 

1.84 

21.5 

1340 

1030 

23.1 

Heat supply 
inadequate 

15 

Copilot * s 

46 

^64 

3.3 

4.95 

27.1 

1280 

935 

27.0 

Heat supply 
adequate 

l6 

Copilot ’ s 

47 

63 

1 

1.0 

L 

1.03 

23.3 

1120 

920 

17'. 8 

Heau supply 
: ad,6jquate 

IT 

Pilot’s 

4 l 

1 

, 51 . I 

r 1 

i 0 

i 

0 

20.5 

730 

430 

41.0 

Heat supply 
adequate 

18 

Pilot’s 

38 

50 

0 

0 

30.0 

860 

630 

26.8 

Heat supply 

slightly 

inadeqviate 

19 

Pilot's 

39 

51 

.1 

.07 

28.6 

875 

630 

28.0 

Heat supply 

slightly 

inadequate 

20 

Copilot’s 

39 

53 

1.1 

1.43 

30.1 

840 

670 

20.0 

Heat supply 
adequate 

21 

Copilot ’ s 

39 

53 

2.2 

1.78 

29.5 

830 

670 

19.3 

Heat supply 
adequate 

22 

Copilot ’ 3 

39 

■55 

3.6 

3.45 

30.9 

885 . 

680 

23.2 

Heat supply 
adequate 

23 

Copilot ’ s 

38 

54 

4.9 

3.12 

28.8 

855 

680 

20.5 

Heat supply 
adequate 

24 

Copilot ’ s 

35 

50 

0 

0 

31.4 

J 

570 

490 

14.0 

Heat supply 

slifShtly' 

inadeqviate 


h* Included the effect of ed.ge losses calculated by the emperical equation 
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TABLE. IS.- A CALCULATION OF TES AIIOUNT CF HATER 
BITEECEPTED BT A FL/iT-PLATE PATEL UCPNG TEE 
DROP-SIZE DISTRIBUTION. ICING COTEITTCN 12 


Liquid 

water 

content 

(percent) 

Liquid L. . . 

T Distri— 

water u , . 

, . but 1 on- 

content i T? 

fm/m3) I ^ 

an 10*^ 
(cm) 

(percent) 

w 

(lb /in*) 

5 

0.07 

1 

j 0.23 

2.19 

0 

0 

10 

.14 

1 

4.18 

0 

0 

20 

.28 

I .65 

6.18 

0 

0 

30 

.42 

1 1.00 

P.5 

2 

0.571 

20 

.28 

1 1.48 
1 

l4.08 

. 13 

2.47 

10 

.14 

1 2.00 

1 

19.0 

30 

2.855 

5 

.07 

■ 2.71 

29.8 


2.24 

100 

1 

0 

1 — 1 

1 

1 


1 

8.136 . 


cp = 2.66 10 

Based on mean— effecti-ve drop diameter of 19 microns 
ao = 0.00099 om 

w = 1.91 It /hr ' : 

q - 890 Btu/hr, sq ft 

Based on drop— size distribution E 
V = 8,lL Ih/hr 

q = 1035 Btu/lir, sq ft 

Experimental data 
q = 1235 Btu/hr, sq ft 

Water drop-size distribution and percentages of 

liquid water content were obtained from reference 4. 
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TABLE X.- CALCULATED MIHIMQM VJIIOSfflELD HEAT REQUIEEMEIWS FOR A 
MODERATE ICING CONDITION IN TWO TYPES OP ICING CLOUDS, 
BASED ON AN ASSUMED DEOP-SIZE DISTRIBUTION E 


Co 


Condi-- 
tion 
mjaaber ' 

1 

Greneral 
cloud 
1 type 

r — — • ] 

' True i 

i 

temper— 1 
ature 
! (°F) 

i 

1 

1 Liquid 

1 vater 
content 
1 (gm/m^) 

i j 

Mo aDr— ! 
effective ; 
drop size | 
(microns) j 

1 i 

True 

airspeed 

feh) 

j Flat— plate 
* windshield 

; <1 
! (Btu/hr, 

i sq ft) 

- L 

Flush ' 

windshield ' 

q ' 

(Btu/hr, 
sq ft) 

V 'type 

' windshield ; 

(Btu/hr, 
sq ft) 

I 

Stratus 

15 

0.5 ] 

' 15 j 

150 

540 

0 

675 ■ 

I 

Stratus 


i 

! 15 

( 

300 

560 

410 

6 

0 

II 

Cumulus 

1 

i ^ 

1 — 

1.0 1 

1 

20 1 

f- - -I- 

150 

J io 4 o 

CO 

0 

0 

1440 

II 

Ciuiiulus 

1 

I ° 

1.0 

1 

£0 ! 
i 1 

O 

o 

1 1630 

1120 

2100 


Pressure altitude, 10,000 ft 
External surface temperature, 32° F 
Panel angle from the fuselage, U5° 
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Figure 1.- The C-46 airplane as equipped for icing research flights in the winter 1945-46 
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electrically heated flush 

WlNIDSHiaO — 


THIN METAL 
HEATING STRIPS 


GUN camera 


NACA 
A-I0854 
I M2 46 


Figure 3.- Adjustable flat-plate windshield and flush windshield panel 
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Figure 4.- Electrically heated flush windshield 
viewed from the exterior of the C-46 aiirplane 


panel, 

cockpit. 






NACA TN No, 1434 


61 



(a) side view. 




(b) Top view. 

Figure 5.- V-type electrically heated windshield 
installed on the C-46 airplane. 
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CO- PlUDTlS HALF V'TYPE va#INOSHIELO%'A^inter \<i<vA'47 



RGUeC 6 - SKETCH SHOv^INGTBE THER.tAOCOUPLE UOCJFHOr^S ON 

VAetOUS tW(N03HlEL0 CONFIGURATIONS USGOOM TUE C-<V€> AIRPLANE 
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Figur© 7#" Singl© psiiiGl f*roni 'tlic V*"*typ© windshiGld installation* Heating 
was provided by transparent, electrical— conduct ing film under outer 

glass layer. 
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Figure 0.- Details of external discharge heated-air ice-prevention system for 
flush windshield of C-46 airplane. 
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VELOCITY OVER , fAP, VA 



TRUE AIB 5 PEED, (V\.p.H 

F10UCE9.- COtv\PAR1SON OF AIRPLANE UELOCITY AND THE iJEUOC\TY 
OVEE FLAT-PLHTE PANEL AND FLUSH -TYPE V^INDSHIElD, 

C-A& airplane 5 WINTER 1945-46 
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reuE speed mPH 

F\6UeE 10.- VAE-IATlON OF O^Y AtP HEAT TeAfSSFER. COEFF\ClElST WITH A\R. SPEED 
AMO altitude on FLAT PLATE PANEL SET 30® F60M THE TANGENT 
TO THE FUSELAGE. C-A6 A IBPLANE ^ W INTE6 
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RGueETW.- VAe\AT\ON Or D<SY Aisi PiEAT TRANSrce. COErRG\ENT WITHAVB. 
SPEED AtSO ALTITUDE ON THE FLAT-PLATE PANEL SET AS “ reOiN\ THE 
tangent to TWE fuselage. C-A<o plane , UJINTEB \M5-4S 


-<1 
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TC.ue fMR, 5 PEED 

riGUGElE - VAR.\AT10N OF OBV A\B HEAT TBANSFeC COEFFICIENT WITH A\B SPEED 
ANO ALTITUDE ON THE FLAT-PLATE RC\NEL SET 60“ FBOI^ THE TANGENT 
TO THE fuselage. C-AG PLANE , vAJ INTER 1‘5A5'-A6, 


DO 
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fa) In icing conditions. Time, 1549. 



(b) Immediately after leaving icing conditions. Time, 1551. 

Figure 13.- Flush windshield of C-46 airplane during icing 
condition 6, table I, with inadequate heat supply to the 
external discharge system for the pilot’s panel. 
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TIME 

PRESSURE ALTITUDE 

TRUE AIR SPEED-— 

AMBIENT AIR TEMPERATURE — 

LIQUID WATER CONTENT 

MEAN EFFECTIVE DROP SIZE- 

HEAT SUPPLIED TO PANEL 

HEAT SUPPLIED TO WINDSHIELD 


1650 TO 1720 

10,300 TO 10,600 FEET 

170 TO 150 MPH 

20 TO I7.5*F 

.15 TO .3 GRAM/M^ 

50 TO 20 MICRONS 

1205 BTU/HR FT^ 

525 BTU/HR FT^ 


Figure 14.- Sketch based on photographs and flight engineer's notes 
of ice accretions on nose and windshields of C-46 airplane after 
icing condition 9, table I. Panel angle 60° with the fuselage. 
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(a) In icing conditions, start of removal. Time, 1703. 



(b) In i cing conditions, removal continued. Time, 1706. 

Figure 15.- Windshield of C-46 airplane during icing condition 9, 
table I, showing removal of ice accretions from the electrically 
heated flush panel, right side, and the external discharge area, 
left side. 
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(c) In clear air, removal completed. Time, 1720. 



(d) 


Ice formation on inch- diameter rod during the 


icing run. 


Figure 15. - Concluded. 
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Figure 16.- Windshield of C-46 airplane after a flight in 
icing conditions in which there was no tendency for ice 
to form on the windshield. Ice accretions shown formed 

on --inch-diameter rods beside windshields; icing 
8 

conditions 10 and 11, table I# 
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FLIGHT NUMBER 31 
PRESSURE ALTITUDE 
TRUE AIR SPEED 
AMBIENT AIR TEMPERATURE 
LIQUID WATER CONTENT 
MEAN DROP SIZE 
HEAT SUPPLIED TO PANEL 
ICING CONDITION NUMBER 


TIME 1116 
8,000 FEET 
160 MPH 
14 •F 
.2 GRAM/M3 
10 MICRONS 
1520 BTU/HR FT^ 




FLIGHT NUMBER 31 
PRESSURE ALTITUDE 
TRUE AIR SPEED 
AMBIENT AIR TEMPERATURE 
LIQUID WATER CONTENT 
MEAN DROP SIZE 
HEAT SUPPLIED TO PANEL 
ICING CONDITION NUMBER 


TIME 1628 
11,650 FEET 
180 MPH 
17 «F 
.5 GRAM/M3 
18 MICRONS 
1240 BTU/HR FT^ 




FLIGHT NUMBER 33 
PRESSURE ALTITUDE 
TRUE AIR SPEED 
AMBIENT AIR TEMPERATURE 
LIQUID WATER CONTENT 
MEAN DROP SIZE 
HEAT SUPPLIED TO PANEL 
ICING CONDITION NUMBER 


TIME 1147 
8,300 FEET 
170 MPH 
18 ®F 
.9 GRAM/M3 
10 MICRONS 
1230 BTU/HR FT2 
3 N A C A 

A-i 0 Q 8 I 


Figure 17.^ Ice accretions formed on the flat-plate windshield of the 
C-46 airplane during icing conditions 1, 2, and 3, table !• Panel 
angle 30^ with the fuselage. 
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FLIGHT NUMBER 2 3 
PRESSURE ALTITUDE 
TRUE AIR SPEED 
AMBIENT AIR TEMPERATURE 
LIQUID WATER CONTENT .6 

MAXIMUM DROP SIZE II 

HEAT SUPPLIED TO PANEL 
ICING CONDITION NUMBER 4 


TIME 1454 
10,350 FEET 
135 MPH 

TO .8 GRAM/M3 
TO — MICRONS 
910 BTU/HR FT* 




FLIGHT NUMBER 23 TIME 1516 

PRESSURE ALTITUDE 11,500 FEET 

TRUE AIR SPEED 145 MPH 

AMBIENT AIR TEMPERATURE -0.5 *F 

LIQUID WATER CONTENT .3 TO 1.2 GRAM/M3 

MAXIMUM DROP SIZE 20 TO — MICRONS 

HEAT SUPPLIED TO PANEL 16 30 BTU/HR FT* 
ICING CONDITION NUMBER 5 




FLIGHT NUMBER 23 TIME 1548 

PRESSURE ALTITUDE 12,700 FEET 

TRUE AIR SPEED 140 MPH 

AMBIENT AIR TEMPERATURE -2 *F 

LIQUID WATER CONTENT .6 TO 1.4 GRAM/M3 

MEAN DROP SIZE 18 TO 44 MICRONS 

HEAT SUPPLIED TO PANEL 172 0 BTU/HR FT* 
ICING CONDITION NUMBER 6 naca 

A 10 982 

I 2-1 1-46 


Figure 18.- Ice accretions formed on the fiat-plate windshield of the 
C-46 airplane during icing conditions 4, 5, and 6, table I. Panel 
angle 45^ with the fuselage. 
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Figure 20.- V-type windshield of C-46 airplane during flight 
in icing conditions showing unsuccessful ice-removal 
operation of exteinal discharge jets over pilot's panel, 
with heat supply of 18,000 Btu per hour. 
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Figure 21.- Typical ice accretion on V-type windshield 
without any heat supplied. 



Figure 22.- Ice accretions formed on V-type windshield of 
the C-46 airplane during icing conditions 15 and 16. 
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(b) In clear air, ice removed# Time, 1215* 

Figure 23.- Ice accretions formed on V-type windshield of the 
C-46 airplane during icing condition 17, 




94 NACA TN No. 1434 



reuE Aie. speed f^p.n 

r\GU6E 24 - UOtNDSHlELD SueFACE Te:^4pg;CiArueE ABOVE 
mB'ENT Alfe TEV\P0E.ATUBES IN CLEAE A\E VOITH 
an aveeage bate of weated avr. suppueo of 

I3>500 Btw/hr C-AG AlEPUANE yJlNTEB I'^AS-AG - 
FLUSH TYPE INOAMIElO- 


NACA TN No. 1434 95 



Fiouee 25 -TF^^PE^^JQmJl^E GEADIENTS in the heated aoUNOAEY 
LAYEE.OVE?. THE FLUSH TYPE WINDSHIELD OF TH£ C.-46 AIR.PLANE. 
Pl5ESSUe.E ALTITUDE 12 >000 FEET , TRUE AlR JSPEEO IS'F lAffH. 
HEATED Aie SUPPUED 13,000 BTU/Hie WINTER ISAfiT-Afe 


LOCAL VELOCITY INS^DC 
BOUNOAe.V LAYER. W.,<^/5EC, 


CD 

Oi 


t.o 


.9 


»6 



A 


200 


ISO 


lOO 


s 


y. 




VIA 




4 




© 




Uo W 

O — 0 -252 ZZS 

+ + -250 3 J 5 

A A -£.65 470 

□ — 0 - £50 470 


NATIONAL ADVISORY 
COMMITTEE FOR AERONAUTICS 


50 1 .. — I ♦ — ' » ' I ♦ ' * » ’* — ^ ~ 

O .1 ^ .3 4 .5 .e .7 S .9 1.0 \.25 1-5 

y, DISTANCE FCOKI SUfeFACJEj INCHES 

FIGUI2.E26.- V'ELOCITY PR.OFILES IN THE HEATTO At R. JET 

OVEe THE FLUSH TYPE VJIND3H\ELD OF THE C-4G AlCPLANE- 
PRESSUfeE A^LTTTUOE IO,600FEET. POIMT OF N\EA&OR,ErAEN T 
9 INCHED FeON\ DlSChP^aSE SLOT- lA^mTEfe 1^45-<Lfi 


2.0 


Q 

4 P 00 

7,500 

16,000 

l ©,000 


NACA TN No. 1434 



60 100 \20 \aO 160 180 200 220 

tbue Aie speed 

R6UeE27.- Cor^PAP.lSOIS OF EXPER.WETSTAL ANO CALCULATED VALUES OF DP.V A'S 
HEAT TPATSSFEe CXOEPnClErNT FOE V^ABlOUS ALTITUDES AND STAGNATION POINTS- 
FLAT-PL/ATE: PANETU set so® Feot^TKC TANSEHT to the fuselage. 

C-4<S AlCPLANE , VAJlNTEfe iS-AT-Afi. 


CD 

-C] 


NACA TN No. 1434 


CD 

00 



FlGUeETZe— COtv\PP»^\S>ON OF exPEe\I^EtSTP\U ANO CAUCUUOrTEO VALUES OF DR.V A\lg 
l-\EAT TBANSrce COEmCieUT ATUAG\OV3S AUTITUOES ftNO STP^GNATtON 
POINTS. FLAT- PLATE PANEl. SET 45' FeOfA TAE TAN GENT TO T WE FUSELAGE. 
C-4G AIRPLANE, UUINTEe l'S4«fr-4G 


> 


NACA TN No. 1434 


MFAT TeANSFETe COemC\ENT, h BTU/we. -3QFT. V 



80 lOO \20 \40 leo 180 £00 220 


TRUE Aie SPEED lAP.H 

FIGURE 2a- COmPAe^SON OF EyPE^^V^^ENTAU AMD CALCULATED VALUES OF OEY AVE 
HEAT transfer COEFFVCI EMT FOR UAR\OUS ALTITUDES AND STAGNftT'ON 
POIMTS. flat-plate PANEL SET CO“ rROt^ TWE TATS GENT TO THE FUSELAGE. 

C-4G AIRPLANE , WINTEQ 1945 -AG 

CD 

CD 


NACA TN No. 1434 



FlGUe.E SG- FLOW LINES AROUND PROviECThD AC.EA Of FI-AT-PLAT& RPiNEL 
USEO IN COr^PUTlNG BATE' OF irAPlNGEI'AE'NT- 





100 NACA TN No. 1434 


STA.B 




FlGUeC31.-FLOW; UNES P^aoutso A COI-^PA^^EO WITH Twe PKOjeCTED 

ABEA OF THE C-4fi AlBPLArNE FUSELAGE AND WINDSHIELD > SHOW I MG 
USED IN CO^^PUTINO RATE OF IfAPlNGENNEMT. 


NACA TN No. 1434 101 






i 


‘ W 


! / 

i *• 

• 


) i 

, ", I 


1 





t 





»>,¥ 


■‘IR 



w 


A 


1 


\ 



#'i 





f 


1 





